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1.  INTRODUCTION 


The  objective  of  this  work  is  to  calculate  the  differential  cross  section  for  single  photon  produc¬ 
tion  in  the  scattering  of  an  unpolarized  electron  by  the  Coulomb  potential  of  a  nucleus,  independent 
of  the  final  momentum  and  spin  state  of  the  electron  and  polarization  state  of  the  emitted  photon. 
This  cross  section  was  calculated  some  time  ago  by  R.  L.  Gluckstern  and  M.  H.  Hull*  and  is  referred 
to  here  as  the  Gluckstern-Hull  formula.  Because  of  the  inherent  difficulty  and  tediousness  of  the 
calculation  and  complicatedness  of  the  result  and  because  only  a  brief  outline  of  the  calculation  was 
published,1  this  work  is  commonly  only  referenced  in  the  literature.  This  differential  cross  section  is 
very  essential  to  any  fully  relativistic  and  quantum  mechanical  calculation  of  bremsstrahlung  pro¬ 
duction  as  a  function  of  angle  in  anisotropic  relativistic  beam-plasma  systems.  Therefore,  it  is  impor¬ 
tant  that  the  Gluckstern-Hull  formula  be  both  fully  understood  and  verified. 

Specifically,  the  Gluckstern-Hull  formula  is  given  by1 


d3o/d2Q  dco  =  (Z2e6/8n)(|p/|/|p|co){8m2  sin2  0  (2E2  4  m2)/(|p|V) 

-  2(5E2  4  2EE7  T  3m2)/(|p|2x2)  -  2(|p|2  -  co2)/(|T|  V) 

+  4E//(|p|2x)  4  (L/|p||p/|)[4Em2  sin2  0  (3com2  -  |p|2E')/(|p|  V) 

4  [4E2(E2  4  E72)  -  2m2(7E2  -  3EE7  4  E72)  4  2m4]/(|p|V) 

4  2co(E2  +  EE7  -  m2)/(|p|2x)]  4  (£1/|p7||f|)[4m2/x2 

-  6co/x:  -  2co(|p|2  -  co2)/(|T|2x)]  -  4£2/(|p7|x:)}.  (1) 


Here,  d3o/d2Q  dco  is  the  differential  cross  section  in  electron-nucleus  bremsstrahlung  for  the  emission 
of  a  photon  of  frequency  co/2tt  into  the  solid  angle  d2Q  =  sin2  0  d0  d<f>.  The  angle  0  is  that  between  the 
photon  and  the  incoming  electron  directions,  and  §  is  the  photon  azimuthal  angle.  The  charge  of  the 
nucleus  is  -Ze  where  e  is  the  charge  of  the  electron.  The  electron  rest  mass  is  m.  The  magnitudes  of  the 
initial  and  final  momenta  of  the  electron  are  |p|  and  |p7|,  respectively.  The  initial  and  final  electron 
energies  are  E  and  E7,  respectively.  The  quantities  x,  T,  L,  £,,  and  £2  are  defined  by 


X  =  E  -  |p|cos  0, 

(2) 

T  =  p- 

(3) 

KEE'  -  m2  +  |p||p/|)/(EE/  -  m2  -  £||p'|)]. 

(4) 

£i  =  in  [(|T|  +  |p'|)/(|T|  -  Ip'DL 

(5) 

e2  =  ln[<E'  +  |p/|)/(E/  -  Ip'I)]. 

(6) 

*R.  L.  Gluckstern  and  M.  H  Hull.  Physical  Review.  90  (1953).  1030. 
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where  £  is  the  photon  momentum  (|£|  =  co/c).  Relativistic  units  are  used  throughout  with 
=  c  =  1  and  e2  =  1/137. 

In  section  2,  the  more  well-known  Bethe-Heitler-Sauter  formula2  is  derived  for  the  differential 
cross  section  for  bremsstrahlung  production  in  which  a  photon  of  a  given  frequency  is  emitted  in  a 
given  direction  and  the  secondary  electron  travels  in  a  prescribed  direction,  disregarding  polarization 
effects,  namely. 


d5o/d2Q  d2Q'  dco  =  (Z2e6m2/4rr2)(|p/|/|p||q|4co) 

*  {(|q|2/ xx/m2)(2E2  +  2E'2  -  |q|2)  +  |q|2(l/x  -  1/x')2 

-  4(E/x'  -  E'/x)2  +  (2co|q|2/ m2)(l/x'  -  1/x) 

-  (2co2/m2)(x'/x  +  x/x')},  (7) 


where 


x  —  E  -  (£/|£|)  •  p  =  E  -  |p|  cos  0, 

(8) 

k'  =  E'-  (£/|£|)  •  p'  =  E'  -  |p'|  cos  0', 

(9) 

q=;'  +  iu;. 

(10) 

Here,  d2Q'  =  sin2  O'  d9 '  d<f»',  where  0'  is  the  angle  between  the  scattered  electron  and  emitted  photon 
directions,  and  <f>'  is  the  angle  between  the  component  of  incident  electron  momentum  transverse  to 
the  photon  momentum  and  the  component  of  scattered  electron  momentum  transverse  to  the  photon 
momentum.  In  section  3,  the  integration  of  the  cross  section  equation  (7)  over  the  direction  of  the 
secondary  electron  is  performed,  and  the  Gluckstern-Hull  formula  equation  (1)  is  thereby  verified. 

2.  BETHE-HEITLER-SAUTER  FORMULA 

The  general  amplitude  for  electron-nucleus  bremsstrahlung  with  one  photon  emission  is 
diagrammatically  depicted  in  figure  1.  In  this  process,  an  incident  electron  is  scattered  by  a  nucleus 
of  charge  -Ze  and,  as  a  result  of  the  interaction,  a  single  photon  is  emitted.  In  the  figure,  the  light 
lines,  heavy  lines,  wavy  line,  and  blob  represent  the  initial  and  final  electron,  initial  and  final 
nucleus,  emitted  photon,  and  five-point  interaction  vertex  for  the  process,  respectively.  This  interac¬ 
tion  vertex  is  exceedingly  complicated.  However,  to  lowest  order  in  the  electromagnetic  coupling 
constant,  neglecting  nuclear  recoil  and  treating  the  nucleus  as  an  external  Coulomb  field  source,  the 
amplitude  may  be  approximated  by  the  sum  of  two  more  elementary  amplitudes  as  shown  in 
figure  2.  Here,  the  Coulomb  field  of  the  static  nucleus  is  represented  by  a  dashed  line,  and  the  lowest 


2V.  B.  Berestetskii,  E.  M.  Lifshitz,  and  L.  P.  Pitaevskii,  Relativistic  Quantum  Theory,  Part  I,  Addison-Wesley  Publishing  Co.,  Inc.,  Reading, 

MA  (1971). 
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order  electromagnetic  vertex  is  represented  by  a  dot.  The  two  parts  of  the  approximate  total 
amplitude  correspond  to  the  photon  being  emitted  either  before  or  after  the  Coulomb  interaction. 
The  conditions  for  validity  of  the  Bom  approximation  are  that  both  the  initial  and  final  velocity  (v 
and  vO  of  the  electron  be  sufficiently  big  and  that  the  nuclear  charge,  -Ze,  be  sufficiently  small, 
namely, 


Ze2/v  «  1, 

(11) 

Ze2/ v' «  1, 

(12) 

and 


Z/137  «  1. 


(13) 


Figure  1.  General  amplitude  for  electron-nucleus  bremsstrahlung  with  one  photon  in  final  state. 


Figure  2.  Approximate  amplitude  in  Born  approximation  (lowest  order  in  electromagnetic  coupling  constant) 
neglecting  nuclear  recoil  and  treating  nucleus  as  source  of  external  Coulomb  field 
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By  elementary  S-matrix  considerations2  the  quantum  mechanical  probability  d3JVW  for  the  oc¬ 
currence  of  a  process  in  which  a  single  particle  of  energy  E  is  scattered  in  a  constant  field  in  an  interac¬ 
tion  volume  V  and  produces  any  number  N  of  particles  with  energies  E'  and  momenta  pi  in  the  final 
state  is  given  by2 


d3NW  =  2tt<5(E/  -  E)|M/,|2(2EV)-1  JT  d^'Qn^GE)"1. 

|o*  1 


(14) 


Here,  E,  is  the  total  energy  in  the  final  state,  and  is  the  scattering  amplitude.  The  latter  is  normally 
determined  by  the  Feynman  rules.  The  cross  section  for  the  process  is  given  by 


d3Na  =  d3nV/j, 


(15) 


where  j  is  the  incoming  particle  flux  density.  The  flux 


density  j  is  given  by 


j  =  v/V, 


(16) 


where  v  is  the  velocity  of  the  incoming  particle.  Clearly, 


v  =  |p|/E. 


(17) 


Substituting  equations  (14),  (16),  and  (17)  into  equation  (15),  then 


d3*a  =  2ttc5(E/ -  E)|M/i|2(2|p|)-1  n  d3^(27r)-3(2E')- 


(18) 


For  the  case  of  bremsstrahlung  with  an  electron  and  a  photon  in  the  final  state,  equation  (18) 
becomes 


d6o  =  2tt<5(E/  +  co  -  E)|M/i|2(2|p|)-1(2n)-6(2co)-,(2E/)-1  d3p'  d3£  (19) 


or 


d6o  =  |M,<|:i(8|p|E/co)-M(E/  +  co  -  E)(2n)-5  d3p'  d3£. 


(20) 


V.  B.  Berestetskii,  E.  M,  Lifshitz,  and  L.  P.  Pitaevsku,  Relativistic  Quantum  Theory,  Part  I,  Addison-Wesley  Publishing  Co.,  Inc.,  Reading 
Mj 4  (1971). 
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Since  kM  =  (co,  k)  and  k2  =  0,  then 


l&l 


=  co. 


(21) 


Therefore, 


d3Jt  =  |£|2  d|ic|  d2Q*  =  co2  do  d2Q*,  (22) 


where  d2Q*  is  the  photon  differential  solid  angle.  Also  since  pi  =  (E',  p')  and  p/2  =  m2,  then 

|p'|2  =  E/2  -  m2  (23) 

and 

d3p'  =  Ip'I2  d|p'|  d2Qp/  =  (E/2  -  m2)  d(E'2  -  m2)1'2  d2Qp/.  (24) 

But 

d(E/2  -  m2)1/2  =  (E/2  -  m2)'1/2  E'  dE'.  (25) 

Therefore, 

d3p'  =  (E/2  -  m2)1'2  E'  dE'  d2Qp/  (26) 

or 

d3p'  =  |p/| E/  dE'  d2Qp/.  (27) 

Substituting  equations  (22)  and  (27)  into  equation  (20),  then 

d6o  =  |M/,|2(8|p|E/o))-M(E/  +  co  -  E)(2n)-5|p'|E'  dE'  d2Qp/co2  dco  d2Q*.  (28) 
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Integrating  over  E',  then 


d5o  =  /|M/f|2(8|p|)“1d(E/  +  o)  -  E)(2n)_5|p/|  dE'  d2Qp/co  dco  d2Q*  (29) 


or 


★ 


d5o  —  |M/1|2(|p/|/8|p|)(2n)"5aj  dco  d2Q*  d2Qp/.  (30) 


The  matrix  element  M/,  is  that  corresponding  to  figure  2  and  is  redrawn  in  figure  3  with  momen¬ 
tum  labels.  Here  p  and  p'  are  the  initial  and  final  four-momenta  of  the  electron,  respectively;  k  is  the 
four-momentum  of  the  emitted  photon;  and  q  is  the  four-vector  momentum  transfer  from  the 
nucleus.  By  energy  and  momentum  conservation,  the  last  is  given  by 


q  =  p'  +  k  -  p.  (31) 


Furthermore,  neglecting  the  change  in  proton  energy  due  to  recoil,  then 


q°  ~  0.  (32) 


M 


fi  “ 


Figure  3.  Feynman  diagrammatic  representation  of  matrix  element  M„. 


The  subset  of  Feynman  rules  to  be  used  in  evaluating  the  matrix  element  M«  in  figure  3  is  as 
follows:2 


1.  Associated  with  each  incoming  electron  line  of  four-momentum  p  is  an  amplitude 
u(p),  where  u(p)  is  the  free  electron  momentum  space  wave  function. 


*  Equation  (30)  corresponds  to  equation  (91.2)  in  V.  B.  Berestetskil,  E  M.  Lifshitz,  and  L.  P  Pitaevskii.  Relativistic  Quantum  Theory.  Part  I 
Addison-Wesley  Publishing  Co.  Inc.,  Reading ,  Mv4  (1971)  The  latter  has  a  typographical  error  omitting  the  factor  of  1/8. 

B  Berestetskii,  E.  M.  Lifshitz,  and  L.  P.  Pitaevskii,  Relativistic  Quantum  Theory,  Part  I,  Addison-Wesley  Publishing  Co.  Inc  Reading 

M A  /1 071)  *  '  *'  *’ 
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tron. 


2.  Each  outing  electron  line  isflssociated  with  the  amplitude  u(p')  of  the  final  elec¬ 


3.  Each  vertex  is  associated  with  a  four-vector  -iyM,  where  yM,  fi  =  0,  1,  2,  3,  are  the 
Dirac  gamma  matrices. 

4.  Associated  with  each  incoming  photon  line  of  four-momentum  k  is  an  amplitude 
\/Tn  em,  where  em  is  the  photon  polarization  four-vector,  and  ^  is  the  same  as  the  index  ^  of  the  vertex 
to  which  it  is  connected.  The  Einstein  sum  convention  for  repeated  indices  is  understood. 

5.  Associated  with  each  outgoing  photon  line  is  an  amplitude  n/"5ti  e*,  where  *  denotes 
complex  conjugation. 

6.  Associated  with  each  internal  electron  line  of  four-momentum  p  and  rest  mass  m  is  a 
factor  i(p  +  m)/(p2  -  m2),  where  p  denotes  yHpH. 

7.  The  four-momenta  of  all  lines  meeting  at  an  interaction  vertex  satisfy  four- 
momentum  conservation. 

8.  An  external  field  line  is  represented  by  a  factor  AH(q),  where  the  latter  is  the  spatial 
Fourier  transform  of  the  electromagnetic  vector  potential  and  q  is  the  four-momentum  transfer  from 
the  external  source. 

9.  Diagrams  that  are  identical  after  removal  of  all  photon  lines  must  have  the  same  sign. 

10.  Bispinor  indices  are  arranged  from  right  to  left,  contrary  to  the  arrows. 

11.  For  n  vertices,  there  is  an  overall  factor  in  of  e"/i. 


Therefore, 


Mfi  =  u(p/)(-iyM)i(p/  +  k  +  m)[(p/  +  k)2  -  m2]"1(-iyv)u(p)(4Tt)1/2E*Av(q)(e2/i) 

+  u(p/)(-iy^)i(p  -  k  +  m)[(p  -  k)2  -  m2]“1("iyv)u(p)(4Tr)1/2E*AM(q)(e2/i)  (33) 


or,  by  interchanging  dummy  indices  ^  and  v  in  the  second  term,  then 


Mfi  =  -e2Av(q)(4Tt)1/2E*u(p/)[yM(f /  +  m)(f/2  -  m2)  1yv  +  yv(f  4*  m)(P  -  m2)_1yM]u(p),‘  (34) 


where 


f  =  p  -  k. 

(35) 

{'  =  p'  +  k. 

(36) 

Using  the  definition  of  equation  (35), 

f2  -  m2  =  (p  -  k)2  -  m2  =  p2  -  2pk  +  k2  -  m2. 


(37) 
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But 


P2  =  m2, 

(38) 

k2  =  0, 

(39) 

and  equation  (37)  becomes 


P  -  m2  =  -2pk  =  -2(Eco  -  |p||£|  cos  0), 


(40) 


where  0  is  the  angle  that  the  emitted  photon  makes  with  the  incoming  electron.  For  the  photon. 


|£|  =  (O. 


(41) 


Equation  (40)  becomes 


where 

x  =  E  -  |p|  cos  0. 
Similarly, 


P  -  m2  =  -2xco, 


F2  -  m2  —  2x/o), 


(42) 


(43) 


(44) 


where 


x/  =  E/  -  |p7|  cos  O',  (45) 

The  external  field  of  the  assumed  static  nucleus  has  a  time  component  only,  given  by  the  static 
Coulomb  potential.  Hence,  neglecting  nuclear  recoil, 
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AM(q)  —  4rrZe|q|_2dM0. 


(46) 


Substituting  equations  (42),  (44),  and  (46)  into  equation  (34),  then 


M/,  =  -e24nZe|q|_2(4r01/2£*u  (p')[yM(f 7  +  m)(2x7a))_1y°  -  y°(f  +  m)(2xco)_1y*1]u(p).  (47) 


Defining  the  matrix  variable  QM  by 


Qm  =  yM( p  +  m)(2x/co)“1y°  -  y°(f  +  m)(2xco)-y, 


(48) 


then  equation  (47)  becomes 


M„  =  -e2(47rZe)|q|"2(4Tr)1/2£*u"(p/)QMu(p). 


(49) 


The  Dirac  matrices  y^  satisfy  the  anticommutation  relations 


{r,  y-}  =  2gr, 


(50) 


where  g1"  is  the  Lorentz  space  time  metric. 


(1  0  0  0\ 
0-1  0  0  \ 
0  0-1  0  I 
0  0  0-1/ 


(51) 


Setting  p  =  v  =  0  in  equation  (50),  then 


y°  =  1, 


(52) 


or  setting  p  =  v  =  k,  where  k  =  1,  2,  3, 


y*  =  -1. 


(53) 
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Using  equations  (49)  and  (52),  then 


M;„  =  -e2(4nZe)|q|-2(4n),/2£„u+(p)Q''V°*u(p/) 

=  -e2(4nZe)  |  q|"2(4n) 1  /2E„ut(p)y0y0Q''*y0*u(p/), 


or 


M),  =  -e2(4nZe)|q|-2(4n),/2£„u  (p)Q‘'U(p/), 


(54) 


where 


QV  =  yOQ^yO* 


(55) 


Therefore,  substituting  equations  (49)  and  (54)  in  the  expression  for  the  cross  section,  equation  (30) 
one  obtains  ' 


d5o  =  zv(47t2)-‘(i;2|a)/i;iiq|4)£;£„ 

•Iu(p/)0,,u(p)](u(p)Q1'u(p/)]  d2Q  d2Q'  da).  (56) 


Defining 


T(p,  p',  k)  =  £|T£v[u  (p/)Ql'u(p)][u  (p)Q1'u(p/)], 


(57) 


equation  (56)  is 


d5o  =  Z2e6(47T2)"1(|p/|co/|p||q|4)T(p,  p* ,  k)  d 2Q  d2Q;  dcu.  (58) 


For  an  unpolarized  incident  electron  and  for  the  case  where  the  final  polarization  states  of  the 
electron  and  photon  are  unspecified,  which  is  the  case  of  interest  here,  one  must  average  over  initial 
electron  spin  states  and  sum  over  final  electron  and  photon  spin  states.  The  cross  section  is  then  given 

by 


d5o  =  Z2e6(4Tr2)-,(|p/|co/|p||q|4)T(p,  p',  k)  d2Q  d 2Qf  dco,  (59) 
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where 


T(p,  p',  k)  =  Vzl  I  ZT(p,  p',  k), 

s  s'  h 


(60) 


and  X,  I ,  I  denote  sums  over  initial  electron  spin  states,  final  electron  spin  states,  and  emitted 
photon  polarization  states,  respectively.  The  factor  of  Vi  arises  in  taking  the  average  over  incident 
electron  spin  states  because  there  are  two  possible  spin  states  for  the  electron. 

Since  in  the  present  calculation  one  is  not  concerned  with  the  state  of  polarization  of  the  emitted 
photon,  then  in  summing  over  photon  polarizations  of  the  final  state,  one  uses  the  relation2 


I  (.  .  .  e*ev  -gM„  .  .  .).  (61) 

h 

Also,  in  summing  over  the  final  electron  spin  states  and  averaging  over  the  initial  electron  states  (un¬ 
polarized),  one  uses  the  following  relation:2 


V2I I  [u  (pOQHu(p)][u  (p)QMu(p')l 

s  s' 

=  TrQM(p  +  m)QH(p/  +  m).  (62) 


Substituting  equation  (57)  into  equation  (60)  and  using  equations  (61)  and  (62),  then 


T(p,  p'f  k)  =  V2I  I  I  eXlu(p') 
s  s'  h 

•  QMu(p)][u  (p)Quu(pO] 

or 


T (p,  p' ,  k)  =  -ViTr  gMVQH(p  +  m)  +  m)  (63) 


or 

T (p,  p',  k)  =  -ViTr  QM(p  +  m)QM(p;  +  m).  (64) 

To  proceed,  one  must  reduce  the  trace  T  given  by  equation  (64).  One  first  notes  according  to 
equations  (48)  and  (55)  that 

*v.  B.  Berestetskii.  E.  M  Lifshitz .  and  L.  P.  Pitaevskii,  Relativistic  Quantum  Theory  Part  1,  Addison-Wesley  Publishing  Co..  Inc.,  Reading. 
MA  (1971). 
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(65) 


QV  =  yOQ-y  =  y®(yO+  fV  +  myO*)y^yO*(2x/W)-. 

-  rV+(f+  +  m)yoyo*(2xw)',. 


However,  it  is  true  that2 


y-  =  yYyO, 


(66) 


Using  equations  (66)  and  (52),  then 


■yOyV+yO  —  y0^ywyO^  —  yV 


(67) 


Using  equations  (66),  (52),  and  (50),  it  follows  that 


y°  = 


(68) 


and 


(69) 


thus. 


In  equation  (65),  the  term  y°y°*f  'VV+  can  be  simplified  by  using  equations  (68),  (52),  and  (67); 


y°y°  f'V  y0’  =  y°z f'VVV, 


(70) 


or 


y°y°*f '  *y“*y’>*  =  y°f/y*'. 


(71) 


Similarly, 


y°y°VV+  =  y°VV  =  yV, 


(72) 


r°/fV7  =  yW2f  y°2  =  y“f  y°, 


(73) 


2  V  B.  Berestetskii,  E.  M  Lifshitz,  and  L.  P.  Pitaevskii,  Relativistic  Quantum  Theory,  Part  I,  Addison- Wesley  Publishing  Co.,  Inc.,  Reading, 
MA  (1971). 
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and 


yy+y0  +  y0+  =  yOyV+yO2  __  yVy0 


Then  substituting  equations  (71)  to  (74)  into  equatio: 


IK 


Q"  =  (y°f/yv  +  my°y”)(2jc/co)'1  -  (y“f  y°  +  my“y°)(2xoj)_1. 


(74) 


(75) 


The  function  T  (p,  p' ,  k)  must  be  evaluated.  Substituting  equations  (48)  and  (75)  into  equation 

(64), 


T(p,  p',  k)  =  -Vi  i T,(p,  p',  k), 

where 

T,(p,  p/ ,  k)  =  Tr  yM'  +  m)(2jc,«)’,y#(p  +  mJy'Xf'  +  m)(2x/co)'1y,'(p/  +  m), 
Tj(p,  p',  k)  =  -TryM(f/  +  m)(2x/cu)‘1y°(p  +  m)yM(f  +  m)(2xco)"1y°(p/  +  m), 
T3(p,  p',  k)  =  -Tr  y°(f  +  m)(2jcoo)*,y„(p  +  +  m)(2jc/co)'1yM(p/  +  m), 

T4(p,  p',  k)  =  Tr  y°(f  +  m)(2xco)"1yM(p  +  m)y"(f  +  m)(2x:co)‘1y°(p/  +  m). 

Using  equations  (66),  (68),  and  (69), 

y°py°  =  y°y^pY  =  y°y»y°  p“  =  rt  pH 

or 

y°py°  =  YoP°  -  Yk p*  =  y*pM  =  p, 


(76) 


(77) 

(78) 

(79) 

(80) 


(81) 


(82) 
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where 


p  ' -  =  (e, 


(83) 


Substituting  equations  (52)  and  (82)  into  equation  (77),  then 


Ti(p,  p7,  k)  —  TryM(f7  4-  m)(2x/co)-1(p  4-  m)(f/  4-  m)(2x/co)-1yM(p/  4-  m). 


(84) 


Using  the  cyclic  invariance  of  the  trace  with  equations  (52)  and  (82),  equation  (80)  can  be  rewritten  as 
follows: 


T4(p,  p7,  k)  =  Tr  y°(p7  +  m)y°(f  4-  m)(2xa))_1y^(p  +  m)yM(f  4-  m)(2xa))_1 


(85) 


or 


T4(p,  p7,  k)  —  Tr(p7  +  m)(f  4  m)(2xco)  *yM( p  4-  m)yM(f  4  m)(2xa))_1 


(86) 


or 

A 

T4(p,  p',  k)  =  Try"(f  +  m)(2xa))‘1(p/  +  m)(f  +  m)(2xa))_1yM(p  +  m)  (87) 

or 

T4(p,  p',  k)  =  Tdp',  p,  -k),  (88) 

where,  in  the  last  step,  equation  (84)  with  the  definitions  of  equations  (35),  (36),  (43),  and  (45)  have 
been  used.  By  using  equation  (88),  then 

T,(p,  p',  k)  =  T„(p',  p,  -k).  (89) 

Next,  using  the  cyclic  invariance  of  the  trace,  equation  (79)  may  be  rewritten 

T3(p,  p',  k)  =  -Tr  (p'  +  m)y°(f  +  m)(2xcu)_1yM  (p  +  m)y°(f/  +  m)(2x/o))"1yH.  (90) 
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It  is  also  true  that2 


Titf.S,  .  .  .  a„)  =  Tr(£ _ 5,5,). 


(91) 


Thus,  using  equation  (91),  equation  (90)  may  be  rewritten 


T3(p,  p',  k)  =  Tr  -  yjf'  +  mJUx'o))-1)^  +  m)y(f  +  m)(2xa))-1y°(p/  +  m).  (92) 


Comparing  equation  (92)  with  equation  (78),  then 


T3(p,  p' ,  k)  =  T2(p,  p',  k). 


(93) 


Substituting  equations  (89)  and  (93)  into  equation  (76),  then 


T(P,  p',  k)  =  -y2[T4(p/,  p,  -k)  +  2T2(p,  p',  k)  +  T„(p,  p',  k)] 


(94) 


or 


T(P,  P‘ .  k)  =  T4s(p,  p',  k)  -  T2(p,  p',  k),  (95) 

where 

T4s(p,  p',  k)  =  -V2[T4(p,  p',  k)  +  T4(p',  p,  -k)],  (96) 

and  T2  and  T*  are  given  by  equations  (78)  and  (87),  respectively. 

The  term  T2  must  first  be  evaluated.  By  substituting  equation  (52)  into  equation  (78)  and  using 
equation  (50),  then 

T2  =  -Try^f'  +  m)(2x/co)"1y°(p  +  m)y°yV(f  +  m)(2?cco)"1y0(p/  +  m)  (97) 


or 


2V.  B.  Berestetskii,  E.  M.  Lifshitz,  and  L.  P.  Pitaevskii,  Relativistic  Quantum  Theory,  Part  I,  Addison- Wesley  Publishing  Co.,  Inc.  Reading 
MA  (1971).  •  * 
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T2  =  -Try^f'  4-  m)(2x/a>)'1(p  +  m)(2g°M  -  y**y°)(f  4*  m)(2xco)_1y°(p/  4*  m) 


(98) 


or 


T2  =  Tr[-2y°(f/  4-  m)(2x/a;)"1(p  4  m)(f  4  m)(2xa))“1 *y°(p/  4  m) 

4  yK(f'  4  m)(2x/a>)“1(p  4  m)yMy°(f  4  m)(2xco)_1y°(p/  4  m)],  (99) 

Using  the  cyclic  invariance  of  the  trace  and  equation  (82),  then  equation  (99)  becomes 

T2  =  Trl^f'  4  m)(2x/o>)_1(p  4  m)(f  4  m)(2xco)_1(p/  4  m) 

4  yJS*  4  m)(2x/co)_1(p  4  m)yM(f  4  m)(2xco)_1(p/  4m)].  (100) 

Using  the  following  relations  obtained  from  equation  (50), 


r%  =  4, 

(101) 

Y*h*  —  -2a, 

(102) 

II 

<-Q 

(103) 

in  equation  (100),  then 

T2  =  Tr[-2(f 7  4  m)(2x/co)"1(p  4  m)(f  4  m)(2xo;)‘1(p/  4  m) 

^  £  A 

4  (4F  •  p  -  2mf '  -  2mp  4  4m2)(2x/co)“1(F  4  m)(2xco)_1(p/  4m)].  (104) 

Performing  the  multiplications  and  recalling  that  the  trace  of  the  product  of  an  odd  number  of  gam¬ 
ma  matrices  is  vanishing,2  equation  (104)  becomes 


2V  B  Berestetskii,  E.  M.  Lifshitz,  and  L.  P.  Pitaevskii,  Relativistic  Quantum  Theory,  Part  I,  Addison-Wesley  Publishing  Co,,  Inc.,  Reading, 

MA  (1971). 
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T2  =  (2x/coxco)“1Tr[-f'  pf  p; 


_  A  A  A  A  .A  ^ 

m2f 'p  -  m2f 'f  -  m2f  /p/ 


Orf  rt  rv  Orf 

-  irrpf  -  m2pp' 

-  m2f  p'  -  m4  +  2f'-p(h'  +  m2) 

A  A  A 

-  mf  /(mf  +mp/)-mp(mf 

A 

+  mp;)  +  2m2(fp/  +  m2)]. 


The  trace  in  equation  (105)  is  evaluated  by  using  the  relations2 


Trab  =  4a  •  b 


and 


Tr  abed  =  4[(a  •  b)(c  •  d)  -  (a  •  c)(b  •  d)  +  (a  •  d)(b-  c)]. 


Thus, 


T2  =  (2xK/co2)"1[-4(f/  •  pf  •  p  '  -  V  •  fp  •  p ' 

+  i'  •  P  'p  ’  f)  -  4m2f'  •  p  -  4m2F  •  f  -  4m2f/  •  p '  -  4m2p  • 

f\j  f\j  r\j 

-  4m2p  •  p '  -  4m2f  •  p '  -  4m4  +  8f'  •  pf  •  p' 

+  8m2f/  •  p  -  4m2f/  •  f  -  4m2f/  •  p'  -  4m2p  •  f* 

-  4m2p  •  p'  +  8m2r  •  p'  +  8m4]. 


"V.  B.  Berestetskil,  E.  M.  Lifshitz,  and  L.  P.  Pitaevskil,  Relativistic  Quantum  Theory,  Part  1,  Addison-Wesley  Publishing  Co  Inc 
MA  (1971).  "  " 


(105) 


(106) 


(107) 


f 


(108) 


Reading, 
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Combining  terms,  then  equation  (108)  becomes 

T2  =  (2xx'co’)-‘[-4(f'  •  p; f  •  p - '  -  F  •  fp  •  p - ' 

f\j  f\j 

+  P  *  p  'p  •  f  -  2F  •  p  f  *  pO  -  4m2(f  *  p 

r\j  r\j  f\j  r\j  r\j 

+  F*p/  +  F’f-F,p  +  p  •p/  +  f,p/ 

+  F-f~  +  F-p'  +  r-  p+  p*p'-2?-p')  +  4m4].  (109) 

Using  the  definition  of  p  ,  equation  (83),  it  follows  that 

P  *  P' '  -  P'  *  p  •  (110) 

Substituting  equations  (35)  and  (36)  into  equation  (109),  performing  the  multiplications  using  equa¬ 
tion  (110),  and  combining  terms,  one  obtains 

T2  =  2(xx/co2)“l[-m2p  *  p  +  m2k  *  p  -  m2p'  •  p  '  -  m2k  •  p  7 

r\t 

-  2m2p  *  p'  -  m2xo)  +  m2k  •  k  -  m4  +  m2k  •  p 

r\j  r\j  r\j  r\j  r\j 

+  (p'  *  P  )2  -  (p'  *  p  k  •  p  0  +  k  •  p  p  *  p  7 

r\j  r\j 

+  (p  •  pO2  -  (k'  -  K)wp  ■  p/  -  p'  ■  p  'p  '  p 

<~\j  r\j  r\j  r\j  ^ 

+  p'  *  p'p  *  k-p  •  pk  *  p7].  (Ill) 


To  reduce  equation  (ill)  from  the  expressions  involving  four-vectors  to  ones  containing  explicit 
space  and  time  components,  one  proceeds  as  follows.  By  defining  the  unit  vector  nin  the  direction  of 
photon  emission. 


n  =  k/|k|  =  k/c o. 


(112) 


and  using  equation  (43),  then 


n-  p  -  E-x. 


(113) 
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Similarly,  using  equation  (45)  , 


n  •  p7  =  E7  -  x' .  (114) 

Using  the  definition  of  p  of  equation  (83),  then 

p  •  p  =  E2  +  |p| 2.  (115) 

But  using  equation  (38), 

|p|2  =  E2  -  m2.  (116) 

Substituting  equation  (116)  in  equation  (115),  then 

p  •  P  =  2E2  -  m2.  (117) 

Similarly, 


k  •  p  =  c»)(E  +  n  •  p)  =  c*)(2E  -  k), 

(118) 

p/  •  p 7  =  E/2  +  |p'|2  =  2E/2  -  m2. 

(119) 

k  •  p 7  =  o)(E7  +  n  •  p7)  =  o)(2E7  -  k'), 

♦  ♦ 

(120) 

p  •  p7  =  EE7  -  p  •  p7. 

(121) 

k  •  k  =  2  co2, 

(122) 

k  ■  p  =  co(E  -  n  •  p)  =  cox. 

(123) 

p/  .  p  =  E7E  +  p'  '  P, 

(124) 
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p'  •  p  k  •  p ' '  =  <u(2E'  -  x')(EE'  +  p  •  p'),  (125) 

k  ■  p  p  •  p '  —  (EE'  +  p  •  p')(coE  +  con  •  p) 

=  (EE'  +  p  •  p')(2E  -  x)o>,  (126) 

p'  ■  P  'p  •  p  =  (E'2  +  |p'|2)(E2  +  |p|2)  =  (2E'2  -  m2)(2E2  -  m2),  (127) 

p'  •  p 'p  •  k  =  (E'2  +  |p'|2)(E  +  n  •  p)o>  =  (2E'2  -  m2)(2E  -  x)co,  (128) 

p  •  P  k  •  p '  =  (E2  +  |p|2)(E'  +  n  •  p')co  =  (2E2  -  m2)(2E'  -  x')w.  (129) 

Also,  for  the  elastic  collision  at  hand,  equations  (31)  and  (32)  imply 

E'  -  E  -  a).  (130) 

Substituting  equations  (117)  to  (130)  into  equation  (111)  and  simplifying,  one  obtains 

T2  =  4(wc'o)2)'1{[m2  +  a)2  +  o)(x'  -  x)]p  •  p'  +  (p  •  p')2  -  m2EE'  +  m2cu2 

-  E2E'2  +  cuEE'2  -  ojE2E'  *  coxE'2  +  wx'E2}.  (131) 

The  quantity  p  •  p;  appearing  in  equation  (131)  may  be  rewritten  in  terms  of  the  momentum 
transfer  qas  follows.  Using  equation  (31),  the  space  component  qof  the  four- vector  q  is  given  by 

q=  y  +  t-p.  (132) 

Dotting  the  vector  q  into  itself,  then 

|q|2  =  |p;  +  p|2  =  |p7|2  +  2p;  •  t  +  |k|2 

-  2p  •  p'  -  2p  •  t  +  |p|2.  (133) 

Substituting  equations  (23),  (112),  (114),  (21),  (113),  and  (116)  into  equation  (133)  and  solving  for 
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P  •  p'  =  Yz[ E2  +  E'2  +  co2  -  2m2  -  2co(E  -  E')  +  2co(x  -  x')  -  |q|2]. 


(134) 


It  follows  then  that 

(p  •  p')2  =  Vi[E4  +  2E2E/2  +  6EW  -  4E2m2  -  4E3co  +  4E2E'co 

+  4(E2  +  E/2  +  co2  -  2m2  +  2coE'  -  |q|2  -  2coE)co(x  -  x') 

-  2E2|q|2  +  E/4  +  2E/2co2  -  4E/2m2  -  4EE/2co  +  4E/3co 

-  2E/2|q|2  +  co4  -  4m2a>2  -  4Ea>3  +  4E/a>3  -  2co2|q|2  +  4m4 
+  8m2a>E  -  8m2coE'  +  4m2|q|2  -  8EE'co2  +  4co2E/2 

+  4Ea>|q|2  -  4coE'|q|2  +  4co2(x  -  x')2  +  |q|4].  (135) 

Next,  substituting  equations  (134)  and  (135)  into  equation  (131),  then 

T2  =  (xx/co2)"’{E4  +  2E2E/2  +  6E2co2  -  4E2m2  -  4E3co  +  4E2E'co 
+  4(E2  +  E/2  +  co2  -  2m2  4-  2COE'  -  |q|2  -  2coE)co(x  -  x') 

-  2E2|q|2  +  E'4  +  2E/2co2  -  4E/2m2  -  4EE/2co  +  4E/3co 

-  2E/2|q|2  +  co4  -  4m2co2  -  4Eco3  +  4E'co3  -  2co2|q|2  +  4m4 

+  8m2coE  -  8m2coE/  +  4m2|q|2  -  8EE'co2  +  4co2E/2  +  4Eco|q|2 

-  4coE'|q|2  +  4co2(x  -  x')2  +  |q|4 

+  2[m2  +  co2  4-  co(x'  -  x)][E2  4-  E/2  4-  co2  -  2m2  -  2co(E  -  E') 

4-  2co(x  -  x')  -  |q|2]  -  4m2EE/  4-  4m2co2  -  4E2E/2  4-  4coEE/2  -  4coE2E/ 

-  4coxE/2  4-  4cox'E2}.  (136) 

Simplifying  equation  (136)  and  substituting  equation  (130),  namely, 

E'  =  E  -  co,  (137) 
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then 


T2  =  (xx7 o>2)-1(-8m2E2  +  8m2Eco  -  4E2|q|2  +  4Eo>|q|2  -  2co2|q|2 

4  2m2|q|2  4  |q|4)  -  2(|q|2/co)(x/-1  -  x"1)  4  4(Ex7“1  +  Fx-1).  (138) 

Next  one  proceeds  by  evaluating  T4(p,  p7,  k)  given  by  equation  (87).  By  using  the  cyclic  in¬ 
variance  of  the  trace,  equation  (87)  may  be  rewritten  as 

T4(p,  p7,  k)  =  Tr(f  4*  m)(2xco)"1(p/  4*  m)(f  4*  m)(2xco)'1yM(p  4*  m )yM.  (139) 


Using  equation  (102)  in  equation  (139),  then 

T4(p,  p',  k)  =  Tr(f  4  m)(2xco)“i(p/  4  m)(f  4  m)(2xco)_1(-2p  4  4m),  (140) 

or  multiplying  terms  and  recalling  that  the  trace  of  a  product  of  an  odd  number  of  gamma  matrices  is 
vanishing,  then 

T4(p,  p7,  k)  =  (4x2co2)“,Tr(-2f  4  4m2f  ?7 

4  4m2ff -2m2fp  4  4m2?7f 

-  2m2p7p  -  2m2f  p  4  4m4).  (141) 

Again  using  cyclic  invariance  of  the  trace  to  combine  the  second  and  fifth  terms,  then 

T4(p,  p7,  k)  =  (4x2o>2)“1Tr(-2f  £ 7  f  p 

4  8m2f  £7  4  4m2f  f  -  2m2f  p 

-  2m2p7p  -  2m2 f  p  4  4m4).  (142) 

Next,  using  equations  (106)  and  (107)  and  combining  terms,  equation  (142)  becomes 

T4(p,  p7,  k)  =  (4x2co2)_1(-16f  •  p7  f  *  p  4  8f2  pf  •  p  4  32m2f  •  p7 

4  16m2P  -  16m2f  •  p  -  8m2p7  *  p  4  16m4).  (143) 
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Substituting  equation  (35)  into  equation  (143),  using  equations  (38)  and  (39),  and  combining  terms, 
then 


T4(p,  p7,  k)  =  (4k2co2)_1[“16(p  •  p7-  px  •  k)(m2  -  k  •  p) 

+  32m2p  •  p7  -  16p  •  kp7  •  p  -  32m2k  *  p7 
-  16m2p  •  k  +  16m4].  (144) 

Again,  multiplying  and  combining  terms,  equation  (144)  becomes 

T«(p,  p',  k)  =  4(x2o>2)'1( m2p  *  ?  -  m2p7  -  k  -  m2p  •  k  -  p7  •  kk  •  p  +  m4).  (145) 

One  notes  using  equation  (83)  that 

p  •  p7  =  EE  7  +  p  *  p 7.  (146) 

Substituting  equations  (120),  (123),  and  (146)  into  equation  (145),  then 

T4(p,.  p7,  k)  =  4(x2co2)_1[m2EE7  +  m2p  *  p'  -  m2o>(2E7  -  x7)  -  m2cox  -  co2(2E7  -  x7)x  +  m4],  (147) 

or 

T4(p,  p7,  k)  =  4(x2co2)’1[in2p  •  p' 

+  m2EE7  -  2m2coE7  -  m2co(x  -  x7)  -  2co2E7x 
+  co2xx7  +  m4].  (148) 

Substituting  equation  (134)  into  equation  (148),  then 

T4(p,  p7,  k)  =  4(x2co2)_1  {( m2/2)[E2  +  E/2  +  co2  -  2m2  -  2co(E  -  E7)  +  2co(x  -  x7)  -  |q|2] 

+  m2EE7  -  2m2coE7  -  m2co(x  -  x7)  -  2co2E7x  +  co2xx7  +  m4}.  (149) 
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Combining  terms,  equation  (149)  becomes 


T4(p,  p7,  k)  =  2m2(E  +  E7)2(x2co2r  +  2m2x-2  -  4m2(E  4  E'Hx2©)-* 

-  8E7x-1  4  4x7x-1  -  2m2|q|2(x2co2)“I.  (150) 

Under  the  interchange  (p,  p7,  k)  (p7,  p,-k),  one  sees  from  equations  (130),  (43),  (45),  and  (132) 
that  (E,  E7,  Co,  x,  x7,  q)  “►  (E7,  E,  -co,  x7,  x,  -q).  Therefore,  by  using  equation  (150), 

T4  (p7,  p,  -k)  =  2m2(E  +  E7)2/(x7co)2  4  2m2/x/2  +  4m2(E  4  E7)/x/2co 

-  8E/x7  4*  4x/x7  -  2m2|q|2/x/2co2.  (151) 

Substituting  equations  (150)  and  (151)  into  equation  (96),  then 

T4s  =  -1/2[2m2(E  +  E7)2/x2o>2  +  2m2(E  4  E7)2/x7V  +  2m2/x2  +  2m2/x/2 

-  4m2(E  4  E7)/x2co  4  4m2(E  4  E7)/x/2co  -  8E7/x  -  8E/x7  4  4x7/x  4  4x/x7 

-  2m2|q|2/x2co2  -  2m2|q|2/x/2co2].  (152) 

Combining  terms  in  equation  (152),  then 

T4s  =  -m2(E  4  E7  -  co)2/ x2co2  -  m2(E  4  E7  4  co)2/x/2co2  4  4(E7/x  4  E/x7) 

-  2(x7/x  4  x/x7)  4  m2|q|2(x-2  4  x7_2)/co2.  (153) 

Next,  using  equation  (130)  in  the  first  two  terms  of  equation  (153),  then 

T4s  =  -4m2E/2/co2x2  -  4m2E2/co2x/2  4  4(E7/x  4  E/x) 

-  2(x7/x  4  x/x7)  4  m2|q|2(x“2  4  x7"2)/co2  (154) 
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or 


T  =  -(4m2/co2)(E7/x  -  E/x7)2  -  (8m2/a>2)(EE7/xx7)  +  4(E7/x  +  E/x7)  -  2 (x7/x  +  x/x7) 

4s 

+  m2|q|2(x-1  -  k'-')2/cj2  +  2m2|q|2/o>2xx2.  (155) 

Returning  to  equation  (95)  and  substituting  equations  (138)  and  (155),  then 

T  =  -(4m2/a>2)(E2/ x  -  E/x2)2  -  (8m2/cq2)(EE2/xx2)  +  4(E2/x  +  E/x2)  -  2(x'/x  +  x/x2) 

+  (m2|q|2/a>2)(x“'  -  x2'1)2  +  2m2|q|2/a>2xx2  -  (l/xx2a>2)(-8m2E2  +  8m2Ea>  -  4E2|q|2 
+  4Ew|q|2  -  2cq2|q|2  +  2m2|q|2  +  |q|4)  +  (2|q|2/cq)(x/-'  -  x")  -  4(E/x2  +  E'/x).  (156) 

One  notes  using  equation  (130)  that 

-4E2|q|2  +  4Ecq|q|2  -  2co2|q|2 

=  (-4E2  +  4 Ecu  -  2cu2)|q|2 

=  [-2E2  -  2(E2  +  co)2  +  4(E2  +  co)co  -  2cu2]|q|2 

=  -(2E2  +  2E/2)|q|2.  (157) 

Substituting  equation  (157)  into  equation  (156),  then 

T  =  (m2/oj2)[(|q|2/xx/m2)(2E2  +  2E/2  -  |q|2)  +  |q|2(x-'  -  x2’1)2 

-  4(E/x2  -  E7x)2  +  (2cq|q|2/m2)(x/-  -  x")  -  (2w2/m2)(x2/x  +  x/x2)].  (158) 

Finally,  substituting  equation  (158)  into  equation  (59),  then 

d5o  =  (Z2ar.2/4Tt2)(|p|2m4/pq4)(dcq/cq)  d2Q  d2Q2  •  [(|q|2/xx2m2)(2E2  +  2E/2  -  |q|2) 

+  |q|2(x-  -  X2")2  -  4(E/x2  -  E2/x)2  +  (2a>|q|2/m2)(x2”1  -  x“) 

-  (2co2/m2)(x2/x  +  x/x2)],  (159) 
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where  r«  is  the  classical  radius  of  the  electron. 


Te  =  e2/m  •  1  =  e2/mc2  =  2.818  x  10  15  meters. 


(160) 


and  a  is  the  fine-structure  constant  (relativistic  units). 


a  =  e2/l  ■  1  =  eVfic  =  1/137.04. 


(161) 


Equation  (159)  agrees  in  form  and  notation  with  equation  (91.7)  of  Berestetskii  et  al.2  It  is  the  well- 
known  Bethe-Heitler-Sauter  formula  for  the  bremsstrahlung  cross  section  in  which  a  photon  of  given 
frequency  is  emitted  in  a  given  direction  and  the  incident  electron  is  scattered  in  a  prescribed  direction 
disregarding  polarization  effects. 


The  cross-section  equation  (159)  can  be  written  in  terms  of  the  scattering  angles  0  and  O'  as 
follows.  One  first  defines 


I  -  (|q|2/xx/m2)(2E2  +  2E/2  -  fq|2)  +  |q|2(xM  -  x'~1)2 

-  4(E /x'  -  E'/x)2  +  (2o)|q|2/m2)(x/“1  -  x"1) 

-  (2co2/m2)(x//x  +  x/xO.  (162) 


Then 


d5o  =  (Z2ar2/4Tr2)(|p/|m4/|p||q|4)(dco/co)  d2Q  d2Q'  X. 


(163) 


Equation  (162)  may  be  rewritten  as 


X  =  Xj  +  X2  +  X3, 


(164) 


3  V.  B.  Berestetskii ,  E  M  lifshitz ,  and  L.  P.  Pitaevskii,  Relativistic  Quantum  Theory.  Part  1  Addison-Wesley  Publishing  Co..  Inc..  Reading. 
MA  (1971). 


30 


where 


I,  =  (m2x;K/)‘1[|q|2(2E2  +  2E/2  -  |q|2)  -  2m2|q|2  +  8m2EE/  +  2co|q|2x: 

-  2co|q|2K/  -  2o)2(?c/2  +  x2)],  (165) 

12  =  Jc-2(|q|2  -  4E/2),  (166) 

13  =  x7-2(|q|2  ~  4E2).  (167) 

The  term  Zj  of  equation  (165)  may  be  written  as 

^1  =  (ir^xx7)-1  |q|2(2E2  +  2E/2  -  | q| 2  -  2m2  +  2c ox  -  2COX7) 

4-  (m2jcx:/)‘l[8m2EE/  -  2co2(x/2  +  x2)].  (168) 

Also,  Z2  of  equation  (166)  can  be  written 

I2  =  (x2m2)_i(|q|2  -  4E/2)m2.  (169) 

But  using  equation  (116), 

m2  =  E2  -  |p|2.  (170) 

Substituting  equation  (170)  into  equation  (169),  then 

I2  =  (ic2m2)_1(|q|2  -  4E/2)(E2  -  |p|2).  (171) 

Substituting 

1  =  sin2  6  +  cos2  6  (172) 


into  equation  (171),  then 


Z2  =  (x2m2)-'(4E/2  -  |q|2)|p|2(sin2  6  +  cos2  d)  +  (7c2m2)-(|q|2  -  4E/2)E2 


(173) 


or 

I2  =  (p2/x2m2)(4E/2  -  |q|2)  sin2  6  +  (l/x2m2)(4E/2  -  |q|2)(|p|2  cos2  0  -  E2).  (174) 

Substituting  equation  (43), 


p|  cos  0  =  E  -  x. 


into  equation  (174),  then 


(175) 


=  (p2/x2m2)(4E/2  -  |q|2)  sin2  0  +  (l/x2m2)(4E/2  -  |q|2)[(E  -  x)1  -  E2]  (176) 


or 


Z2  =  (p2/x2m2)(4E/2  -  |q|2)  sin2  0  +  (l/m2)(4E/2  -  |q|2)  -  (2E/xm2)(4E/2  -  |q|2).  (177) 

Analogously,  Z3  of  equation  (167)  can  be  rewritten  in  the  form 

Z3  =  (p/2/x:/2m2)(4E2  -  |q|2)  sin2  6'  +  (l/m2)(4E2  -  |q|2)  -  (2E//x/m2)(4E2  -  |q|2).  (178) 

Substituting  equations  (168),  (177),  and  (178)  into  equation  (164),  then 

I  =  (l/m2xx/)[|q|2(2E2  +  2E/2  -  |q|2  -  2m2  +  lux  -  lux'  -  2xx'  +  2Ex'  +  2E'x) 

+  8m2EE/  -  2u2(x'2  +  x2)  +  (4E2  +  4E/2)xtc/  -  8EE'V  -  8E'E2x] 

+  (p/2/x/2m2)(4E2  -  |q|2)  sin2  6'  +  (|p|2/?c2m2)(4E/2  -  |q|2)  sin2  0  (179) 

or 


32 


where 


X  =  X/m 2xx;  +  (|p/|2/x/2m2)(4E2  -  |q|2)  sin2  0' 
+  (|p|2/x2m2)(4E/2  -  |q|2)  sin2  0, 


(180) 


» 


*  X  =  Xj  +  |q|2(-2m2  +  2cox  -  2c ox'  -  2xxl  +  2Ex7  +  2E;x)  +  8m2EE/ 

-  2a>2(x/2  +  x2)  +  (4E2  +  4E/2)xx7  -  8EE'V  -  8E/E2x/ 

I.  =  |q|2(2E2  +  2E/2  -  |q|2). 


(181) 

(182) 


Using  equation  (133),  |q|2  can-  be  expressed  in  terms  of  the  scattering  angles  0  and  G/  as  follows.  The 
momentum  vector  p  can  be  written  in  terms  of  components  p^  and  p^  perpendicular  and  parallel  to 
the  photon  momentum  vector  k,  respectively;  thus. 


Clearly, 


♦  _  ♦  .  * 

P  “  P±  +  P[|- 


Pii  =  (p-k/|ic|2)L 


(183) 


(184) 


Substituting  equation  (184)  into  equation  (183),  then 


P  =  Pl  +  (p-  ic/|ic|2)t 


(185) 


and  by  definition 


Analogously, 


PL  •  £  =  0. 


P'  =  Pi  +  ■  t/|C|2)t 


(186) 


(187) 
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and 


pi  •£  =  o. 


(188) 


Also,  denoting  the  angle  between  px  and  hi  by  •)>',  then 


pi  -p1  —  IpiIIpII  cos  <t>/. 


(189) 


Denoting  the  angle  between  the  incident  electron  momentum  p  and  the  emitted  photon  momentum  k 
by  9  and  that  between  the  scattered  electron  momentum  p'  and  £  by  9',  then 


Ipi  1  =  |p|  sin  9, 

(190) 

P-  £  =  |p|  |£|  cos  9, 

(191) 

|pil  =  Ip'I  sin  9', 

(192) 

p'  •  £  =  Ip'I |j<|  cos  9'. 

(193) 

Using  equations  (185)  and  (187),  then 


or 


P  ■  p'  —  [pi  +  (p  •  ic/|ic|2)k] 

■  [pi  +  (p'  ■  C/|i<|2)ic] 


(194) 


p-  p'  =  pi  ■  pi 


■r 


* 
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+  (p-k/|£|2)pi  •  t 
+  (^  •  ic/|ic|2)pj_  •  £ 
+  (p-^)(p'-t)/|C|2. 


(195) 


Substituting  equations  (189),  (186),  (188),  and  (190)  to  (193)  in  equation  (195),  then 


p-  p'  =  |p|  sin  0  |p;  |  sin  9/  cos 

+  (fP||k|/|k|2)  cos  9  Ip7  |  |lc|  cos  9'  (196) 


or 


p •  p/  =  |p|  |^|  (cos  9  cos  0'  +  sin  9  sin  0;  cos  <t>7).  (197) 

Therefore  substituting  equations  (21),  (193),  (191),  and  (197)  into  equation  (133),  one  obtains 

|q|2  =  |p|2  +  |j/|2  +  o)2  -  2|p|co  cos  9  +  2 |p/ |a>  cos  0' 

-  2|p|  I^Kcos  9  cos  9 '  +  sin  9  sin  9(  cos  (198) 

Then  substituting  equation  (198)  into  equation  (182), 

—  [  Ip|2  +  Ip"!2  +  co2  -  2|p|co  cos  0  +  2|p/|co  cos  0 1 

-  2|p||p/|(cos  0  cos  9'  +  sin  0  sin  0'  cos  4>7) ](2E2  +  2E/2  -  |q|2)  (199) 


or 

X,  =  -2|p|  1^1  sin  0  sin  07  cosftfE2  +  2E/2  -  |q|2)  +  (|p|2  +  |pf  +  co2 

-  2|p|co  cos  0  +  2|f/|o>  cos  9'  -  2|p|  |f/|  cos  0  cos  07)(2E2  +  2E/2  -  |q|2).  (200) 

Next  substituting  equation  (200)  into  equation  (181),  then 

I  =  -2|p|  |^|  sin  0  sin  9'  cos  (2E2  +  2E/2  -  |q|2)  +  |q|2  Z2  +  21 3,  (201) 


where 


Z2 - 2m2  +  2c ox  -  2COX2  -  2xx'  +  2Ex2  +  2E2x  -  |p|2  -  Ip7!2  -  co2 

+  2|p|cu>  cos  0  -  21^0)  cos  O'  +  2|p|  Ip7!  cos  0  cos  6',  (202) 

^3  =  (E2  +  E/2)(|p|2  +  Ip7!2  +  co2  -  2|p|co  cos  0  +  2^10;  cos  O' 

-  2|p|  |p'|  cos  0  cos  00  +  4m2EE/  -  co2(x/2  +  x2) 

+  2(E2  +  E/2)xx2  -  4EE/V  -  4E2E2x.  (203) 

Substituting  equations  (43)  and  (45)  in  equation  (202),  then 

Z2  =  -2m2  +  2co(E  -  |p|  cos  0)  -  2co(E/  -  \}/\  cos  00  -  2(E  -  |p|  cos  0)(E'  -  |^|  cos  00 
+  2E(E'  -  |p'|  cos  00  +  2E'(E  -  |p|  Cos  0)  -  |p|2  -  |^|2  -  co2 

+  2|p|cp  cos  0  -  2|p/|co  cos  O'  +  2|p|  Ip7!  cos  0  cos  0'  (204) 


or 


I2  =  -2m2  +  2cjE  -  2co|p|  cos  0  -  2EE'  +  2E\b'\  cos  O'  -  2|p|  Ip7!  cos  0  cos  O' 

-  |p|2  -  2ojE'  +  2cl)\p'\  cos  O'  +  2EE'  +  2E/|p|  cos  0 
+  2|p|  Ip'I  cos  0  cos  0'  -  Ip'l2  +  2co|p|  cos  0 

+  2E2E  -  2E|p/|  cos  O'  -  to2  -  2u\p'\  cos  O'  -  2E'|p|  cos  0.  (205) 

Substituting  equation  (130)  into  the  8th  and  15th  terms  of  equation  (205),  then 

Z2  =  -2m2  +  co2  +  2E2  -  2wE  -  |p|2  -  Ip7!2.  (206) 

Further,  using  equation  (130), 


co  =  E  -  E',  (207) 

and  substituting  equation  (207)  in  the  second  and  fourth  terms  of  equation  (206),  then 
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I2  =  -2m2  +  E2  +  E/2  -  |p|2  -  fpf* 


(208) 


Next,  substituting  equation  (116)  and,  analogously, 

|^|2  =  E/2-m2 


into  equation  (208),  then  finally  • 


(209) 


X2  =  0.  (210) 

Also,  substituting  equations  (43)  and  (45)  into  equation  (203)  and  multiplying  terms,  then 

X3  =  E2|p|2  +  E2|J/|2  +  E2co2  -  2E2|p|co  cos  0 
+  2E2|p/|co  cos  9/  -  2E2|p|  Ip7)  cos  0  cos  9/ 

+  E/2|p|2  +  E/2|^|2  +  E/2co2  -  2E/2|p|co  cos  0 
+  2E/2|p/|co  cos  0'  -  2E/2|p|  1^1  cos  0  cos  0; 

+  4m2EE/  -  co2E/2  +  2co2E/|J>/l  cos  9'  -  co2^2  cos  29' 

-  co2E2  +  2co2E|p|  cos  9  -  co2|p|2  cos  29  +  2E3E7 

-  2E3|p/|  cos  9/  -  2E2E/|p|  cos  9  +  2E2|p|  |{/|  cos  9  cos  9 ' 

+  2EE/3  -  2EE/2|p/|  cos  0'  -  2E/3fP|  cos  9 

+  2E/2|p|  |^|  cos  9  cos  0'  -  4EE/3  +  4EE/2|J/|  cos  0' 

-  4E7E3  +  4E/E2|p|  cos  9.  (211) 

Substituting  the  trigonometric  identities 

cos2  6K  =  1  -  sin2  9/  (212) 


and 
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cos2  0  =  1-  sin2  0 


(213) 


into  the  16th  and  19th  terms  of  equation  (211)  and  combining  terms,  then 

I3  =  co2|p|2  sin2  0  +  co2|p/|2  sin2  61  +  X4,  (214) 

where 

X4  =  E2|p|2  +  EV|2  +  E/2|p|2  +  E/2|p/|2  +  4m2EE/  -  co2|J/|2 

-  co2|p|2  -  2E3E;  -  2EE/3  +  cos  0  (-2E2|p|o>  -  2E/2|p|o> 

+  2co2E|p|  -  2E2E/|p|  -  2E/3|p|  +  4E/E2|p|) 

+  cos  O'  (2E2|p/|o)  +  2E/2|p/|o>  +  2co2E/|p/|  -  2E3|J/| 

-  2EE/2|p/|  +  4EE/2|p/|).  (215) 

Next,  substituting  equation  (207)  into  equation  (215)  and  combining  terms,  then 

Z4  =  4m2EE/  +  2EE/|p/|2  +  2EE'|p|2  -  2E3E'  -  2EE/3.  (216) 

Further,  substituting  equations  (116)  and  (209)  into  equation  (216),  then 

=  0.  (217) 

Then  substituting  equation  (217)  into  equation  (214), 

X3  =  co2|p|2  sin2  0  +  co2|p'|2  sin2  07.  (218) 

Therefore,  substituting  equations  (210)  and  (218)  into  equation  (201),  then 

X  =  -2|p|  Ip7!  sin  0  sin  0'  cos  <t>;  (2E2  +  2E/2  -  |q|2) 

+  2co2|p|2  sin2  0  +  2co2|p/|2  sin2  0'.  (219) 
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Then  substituting  equation  (219)  into  equation  (180), 


I  =  m-2[(|J/|2/x/2)(4E2  -  |q|2)  sin2  O'  +  (|p|2/x2)(4E/2  -  |q|2)  sin2  0 
+  (2a>2/xx7)(|p|2  sin2  0  +  |p/|2  sin2  07) 

-  (2|p|  |p7|/xx7)(2E2  +  2E/2  -  |q|2)  sin  0sin  07  cos  <t>7].  (220) 

Finally,  substituting  equation  (220)  into  equation  (163), 

d5o  =  (Z2ari/4n2)(|p/|m2/|p|  |q|4)  (dco/co)  d2Q  d2Q7 

•  [(|p/|2/x/2)(4E2  -  |q|2)  sin2  O'  +  (|p|2/x2)(4E72  -  |q|2)  sin2  0 
+  (2co2/xx7)(|p|2  sin2  0  +  |f/|2  sin2  07) 

-  (2|p|  |p/|/xx7)(2E2  +  2E/2  -  |q|2)  sin  6  sin  O'  cos  <}>7].  (221) 

The  photon  and  recoil-electron  differential  solid  angles  are  given  by 

d2Q  =  sin  0  dO  d<}>  (222) 

and 


d2Q7  -  sin  0'  dO'  df,  (223) 

respectively.  If  one  integrates  over  the  photon  azimuthal  angle  <{>,  the  corresponding  cross  section  is 
given  by 

d4o  =  /  (dso/d<}>)  d$.  (224) 

o 

Substituting  equation  (221)  into  equation  (224)  and  using  equations  (222)  and  (223),  then 
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d4°  “  (Z2arJ/2n)(|p'|m2/|p|  |q|4)  (dco/co)  sin  0  d0  sin  O'  d6/  d§' 

*  [(|J/|2/k/2)(4E2  -  |q|2)  sin2  O'  +  (|p|2/x:2)(4E/2  -  |q|2)  sin2  0 
+  (2co2/kk/)(\p\2  sin2  0  +  [p7!2  sin2  O') 

-  (2|p| |p'|/joc/)(2E2  +  2E/2  -  |q|2)  sin  0  sin  O'  cos  (225) 

Equation  (225)  agrees  with  equation  (91.8)  of  Berestetskii  et  al/2  wherein  the  notation  for  <f>  and  §'  is 
interchanged  from  that  used  here.  The  differential  cross  section  d5o  of  equation  (221)  is  the  probabili¬ 
ty  (in  Born  approximation)  of  single  photon  emission  into  solid  angle  d2Q  and  electron  recoil  into 
solid  angle  d2Q'  per  unit  electron  flux.  The  functions  x,  k',  and  |q|2  are  given  by  equations  (43),  (45), 
and  (198),  respectively.  The  quantities  E  and  |p|  are  the  energy  and  the  momentum,  respectively,  of 
the  incident  electron;  E'  and  [p7 1  are  the  energy  and  the  momentum,  respectively,  of  the  scattered 
electron;  O'  is  the  polar  angle  made  by  the  scattered  electron  momentum  with  the  photon  momentum; 
o)/2tt  is  the  frequency  of  the  emitted  photon;  0  is  the  polar  angle  between  the  incident  electron 
momentum  and  the  emitted  photon;  <t>7  is  the  angle  between  the  component  of  the  incident  electron 
momentum  transverse  to  the  photon  momentum  and  the  component  of  the  scattered  electron 
momentum  transverse  to  the  photon  momentum;  <j>  is  the  azimuthal  angle  of  the  photon  momentum; 
-Ze  is  the  charge  of  the  nucleus;  cr  is  the  fine  structure  constant  given  by  equation  (161);  and  r  is  the 
classical  radius  of  the  electron  given  by  equation  (160). 


3.  GLUCKSTERN-HULL  FORMULA 


In  this  section,  the  Bethe-Heitler  cross  section,  equation  (221),  is  to  be  integrated  over  all  direc¬ 
tions  of  the  recoil  electron  to  obtain  the  Gluckstern-Hull  formula,  equation  (1).  The  latter  is  the  dif¬ 
ferential  cross  section  in  Born  approximation  for  single  photon  emission  in  the  scattering  of  an  elec¬ 
tron  by  the  Coulomb  potential  of  a  nucleus,  independent  of  the  direction  of  recoil  of  the  final  state 
electron.  Also,  the  calculation  here  is  specialized  to  the  case  of  an  unpolarized  incident  electron,  and 
the  cross  section  is  summed  over  all  polarization  states  of  the  emitted  photon  and  all  final  spin  states 
of  the  recoil  electron. 

Thus,  the  Gluckstern-Hull  cross  section  is  given  by 


d3o  =  / (d5o/d2Q/)  d2Q', 


(226) 


where  the  integral  is  over  all  possible  angles  of  electron  recoil.  By  substituting  equation  (221)  into 
equation  (226),  then 


*V.  B.  Berestetskii,  E,  M.  Lifshitz,  and  L.  P.  Pitaevskii,  Relativistic  Quantum 
MA  (1971). 


Theory,  Part  1,  Addison-Wesley  Publishing  Co.,  Inc.,  Reading, 
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d3o  =  (Z2arf/2n)(|p/|m2/|p|)  (dco/co)  d2Q  I(|p|,  k,  Ip'l), 


(227) 


where 


I(|p|,  k  Ip'l)  =  / (d2Q//2n)|q|_4[(|p'|2/x/2) 

•  (4E2  -  |q|2)  sin2  9'  +  (|p|2/x2)(4E/2  -  |q|2)  sin2  6 

+  (2co2/jcjc/)(|p|2  sin2  0  +  |f/|2  sin2  O')  -  2(|p|  Ip^xx') 

*  (2E2  +  2E/2  -  |q|2)  sin  0  sin  O'  cos  $']. 


(228) 


One  proceeds  to  calculate  I(|p|,  k,  |J/|).  Equation  (228)  is  first  rewritten  in  the  equivalent  form 


I  =  J(d2Q//2n)|q|“«{  |q|2[-(|p'|2/x/2)  sin2  O'  -  (|p|2/x2)  sin2  0 
+  2(|p|  Ip'l /xrx:7)  sin  0  sin  O'  cos  $  ]  +  4(|p'|2E2/x/2)  sin2  0' 

+  4(|p|2E/2/x2)  sin2  0  +  2(w2/xx')(\p\2  sin2  0  +  |pf  sin2  6') 
-  4(|p|  \p'\/xx')(E2  +  E/2)  sin  0  sin  0'  cos  ^  }. 


(229) 


Next,  using  equation  (45), 


cos  0'  =  (E'  -  x')/  Ip'l 


(230) 


Also,  by  trigonometry, 


sin2  O'  =  1  -  cos2  O' . 


(231) 


Substituting  equation  (230)  into  equation  (231),  then 


sin2  0'  =  1  -  (E//|p/|)2  +  2(E/x//fPf )  -  */2/|p12 


(232) 


41 


Also,  from  equation  (198),  it  follows  that 


sin  0  sin  0'  cos  f  =  (Vz  |p|  |p|)(|p|2  +  |pf  +  .  |*ja 

-  2|p|co  cos  0  +  2|p'|co  cos  6'  -  2|p|  |p'|  cos  0  cos  0' 

Next,  substituting  equations  (232)  and  (233)  into  equation  (229),  then 

I  =  /  (d2Q//2n)  { |q|-2  [-|p'|2/ x'2  +  E/2/x/2  -  2E'/x'  +  1 

-  (|p|2/x2)  sin2  0  +  (xx'j-'flpl2  +  |J/|2  +  co2  -  |q|2 

-  2|p|cp  cos  0  +  2cpE2  -  2cpx2  -  2|p|E'  cos  0  +  2|p|x2  cos  0)  ] 

+  |q|'4[4(|p/|2E2/x/2)  -  4(E/2E2/x/2)  +  8(E/E2/x/)  -  4E2 

+  4(|p|2E/2/x2)  sin2  0  +  2(co2|p|2/xx/)  sin2  0  +  2(cp2|J>'|2/xx/) 

-  2(cp2E/2/xx/)  +  4(E'cp2/x)  -  2(coV/x)  -  2(xxO"'(E2 

+  )(|p|2  +  | P'12  +  CP2  -  |q|2  -  2|p|co  cos  0 

+  2cpE'  -  2cpx'  -  2|p|E'  cos  0  +  2|p|x/  cos  0)  ]  }. 

Equivalently, 


1  =  f  (d2Q//2n)  {  |q|~2[  (x/2)_1(-|p/|2  +  E/2)  +  (xx/)',(-2E/x  +  |J>|2 

+  lf/12  +  co2  -  |q|2  -  2|p|cp  cos  0  +  2cpE2  -  2pE2  cos  0)  +  x‘‘(-2cp  +  2| 
+  1  -  (|p|2/x2)  sin2  0]  +  |q|  4 [ x/  2(4|p/|2E2  -  4E/2E2)  +  (x'x)-1  {8E'E2x 
+  2cp2|p|2  sin2  0  +  2cp2|p/|2  -  2cp2E/2  -  2(E2  +  E/2)(|p|2  +  |J/|2  +  cp2 

-  |q|2  -  2|p|cp  cos  0  +  2cpE'  -  2|p|E'  cos  0)}  +  4x‘'(E2  +  E/2)(cp  -  |p|  c< 

-  4E2  +  4(|p|2E/2/x2)  sin2  0  +  4E/cp2/x  -  2cpV/x]  }. 

Using  equation  (209),  one  notes  that 


-|pf  +  E/2  =  m2 


).  (233) 


(234) 


p|  cos  0) 

is  0) 

(235) 

(236) 
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and 


4|p/|2E2  -  4E/2E2  =  -4m2E2. 


(237) 


Substituting  equations  (236)  and  (237)  into  equation  (235),  then 


I  =  f  (d2Q'/2n)  { (m2(|q|2x/2)~‘  +  (|q|2xx')-‘ [-2E'x  +  |p|2  + 

+  w2  -  2|p|ol>  cos  0  +  2cuE'  -  2|p|E'  cos  0  +  2(E2  +  E/2)]  -  (xx')-> 

-  |q|-2 [2x_1(cl>  -  |p|  cos  0)  +  |p|2x“2  sin2  0  -  l]  -  4m2E2|q|~4x2'2 
+  (|q|4xx2)-1  [8E2E2x  +  2cp2|p|2  sin2  0  +  2o)2fP2|2  -  2cu2E/2  -  2(E2  +  E/2)(|p|2 
4-  Ijp2!2  4-  a>2  -  2|p|cj  cos  0  4-  2a)E2  -  2|p|E/  cos  0)]  +  |q|‘4[4x_1(E2  +  E/2) 

•  (a>  -  |p|  cos  0)  -  4E2  +  4x“2|p|2E/2  sin2  0  +  4E/a)2x"1  ]  -  2cj2x/(|q|4x)"1 }.  (238) 


The  terms  to  be  integrated  in  equation  (238)  are  all  of  the  form 


lm,„  =  /( d2QV2n)|q| 


(239) 


where 


{(m,  n)}  =  {(0,  1),  (1,  0),  (1,  1),  (1,  2),  (2,  -1),  (2,  0),  (2,  1),  (2,  2)}. 


(240) 


The  integrals  of  equation  (239)  can  be  cast  in  a  more  transparent  form  as  follows.  Dotting  q  of  equa¬ 
tion  (132)  into  itself,  then 


|q|2  =  |p-^-£|2 
=  |(P-  0  -  pf 

=  |p-  ic|2  +  Ip2!2  -  2f/  ■  (p-  £) 


(241) 


or 
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(242) 


|q|2  =  (|P-£|2  +  fcf) 

•  [l  -  p'  ■  2(p-  lc)  (|p-  lc|2  +  i^i*)-]. 


Defining  vectors  T,  a,  and  b  by 


then  equation  (242)  is 


Also,  using  equation  (45), 


T  =  p  -  £, 

a  =  2T(|T|2  +  |pf)‘*, 
6  =  £(|£|E2)-', 


|q|2  =  (|T2|  +  fcfXi  -  if  ■  D. 


(243) 

(244) 

(245) 


(246) 


—  E'  -  Ipd  cos  0 2 
=  E2(l  -  |p'|E/'1  cos  00. 

Substituting  equation  (193)  into  equation  (247),  then 

x2  =  E2[l-£  •  ^(E'|£|)-‘]. 

Next,  substituting  equation  (245)  into  equation  (248), 

x2  =  E2(l  -?■&). 


(247) 


(248) 


(249) 


Therefore,  substituting  equations  (246)  and  (249)  into  equation  (239), 
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(250) 


I„„  =  /(d2SV2rt)(|T|2  +  fpf  )-"(l  -  P  •  a)- 
•  E'-d  -  p1  ■  by. 


Since  T  is  independent  of  6'  and  ,  it  can  be  taken  outside  the  integral,  and  equation  (250)  becomes 


U  =  (|T|2  +  |;'|2)""E'-"I, 


(251) 


where 


L„  =  J (d2Q//2rt)(l  -  p'  •  a)-“(l  -  p  ■  6)-". 


(252) 


One  proceeds  by  evaluating  the  integrals  of  equation  (252).  Begininning  with  I0i,  one  has 


I01  =  /( d2Q'/2n)(l  -^  •£)-’. 


(253) 


Since  the  photon  momentum  £  has  been  chosen  as  the  polar  axis  in  equation  (193),  then  according  to 
the  definition  of  fi  of  equation  (245),  £  is  along  the  polar  axis.  Therefore, 


p  ■  £  =  Ip2!  |£|  cos  6'. 


(254) 


Expressing  the  element  of  solid  angle  d2Q'  in  equation  (253)  explicitly  in  terms  of  O'  and  <t>'  and 
substituting  equation  (254),  then 


lot  =  /2n(df/27T)/n  sin  0'  dO'  (1  -  fr\\i\  cos  O')'1. 
0  0 


(255) 


One  next  changes  variables  of  integration  by  defining 


x  =  li^l  |£|  cos  0' . 


(256) 


Then 


dx  =  -|p'|  |b|  sin  O'  dO'. 


(257) 
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Then  equation  (255)  becomes  after  integration  over  <j>7 


Ip'I  IbJ 

loi  =  |p'l"1|bl"1  /  dx  (1  -  x)”1.  (258) 

-Ip'llbl 

Performing  the  elementary  integration,  then 

,  .  Ip'llbf 

!oi  =  (M  lbl)~^n  (1  -  x)'1  |  (259) 

Vllbl 


or  finally 


loi  -  (|p'||b|)"1  /n  [(l  +  |p'||b|)(l  -  Ip'llbl)-*)]. 


(260) 


Similarly  evaluating  Ii0,  one  has 


Iio  —  J  (d2Q' /2n)(l  -  f/-  a)-1  =  /  d9'  sin  6'  (1  -  Ip'l  |a|  cos  0')-1, 

0 


(261) 


or  as  in  equations  (255)  to  (260),  then 


Iio  =  (|p/||a|)-1  in[(l  +  |p'|  |a|)(l  -  |p'|a|)-‘]. 


(262) 


Next,  considering  I20, 


I20  =  /  (d2S2V27r)(l  -  ^  •  a)'2 

=  /  d 6'  sin  6'  (1  -  |p'|  |a|  cos  6')~2  (263) 

0 

frill  I 

I20  =  (Ip'l  |a|)"‘  /  ,  ,  dx  (1  -  x)'2 

Vllal 

Ip'llal 

=  (iHfaP-U-*)-*  I.  .  (264) 

-Ip'l  la  I 


or 
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I20  =  2(1-  |a|2  |p/|2)-1. 


(265) 


To  evaluate  I,,,  the  following  identity  is  useful  for  arbitrary  quantities  a  and  ft  independent  of  the  in¬ 
tegration  variable,  namely. 


(ap)  1  =  /  dx  [ax  +  P(1  -  x)]“2.  (266) 

c 

Equation  (266)  is  true  since 

f  dx  [ax  +  p(l  -  x)]“2  =  /  dx  [(a  -  p)x  +  p]~2 
o  o 

=  -(a  -  p)-'[(a  -  p)x  +  PY' 

0 

=  -(a  -  +  (a  -  PY'P-' 

=  (ap)-'.  (267) 

Proceeding  then  with  the  evaluation  of  In, 

1„  =  /  (d2Q//2n)(l-p'  a)-,(l-p' •£)-*.  (268) 

Using  the  identity  equation  (266)  in  the  integrand  of  equation  (268),  then 

Itl  =  /(d2QV2Tr)  /‘dx  [(1  -  p'  *  a)x  +  (1  -  J/  •  b)(l  -  x)]"2,  (269) 


or  interchanging  orders  of  integration  and  rewriting  the  integrand,  then 

In  =  /  dx  /  (d2Q/2rr){l  -  p' •  [ax  +  (1  -  x)b]}“2.  (270) 

0 


Defining  the  vector  g  by 


g  =  xa  +  (1  -  x)b. 


(271) 


then  equation  (270)  may  be  written 
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In  —  /  dx  /  (d2£2/2n)(l  -p'-g)"2. 


(272) 


Clearly,  the  integrals  must  be  invariant  with  respect  to  the  choice  of  the  polar  axis  of  integration.  By 
choosing  the  vector  g  as  the  polar  axis  and  comparing  equation  (272)  with  equations  (263)  to  (265), 
then 


In  —  2/  dx  (1  -  Igplp'l2)"1. 
0 


(273) 


Substituting  equation  (271)  into  equation  (273),  then 


In 


=  2/  dx  [l  -  If/p  |ax  +  6(l-x)|2]- 


(274) 


Equation  (274)  can  be  rewritten  as 


In  =  2/  dx  [l  ■  |^|2|p/|2  -  2x|p'|2(a-  6  -  |6|2 


-x2|p'|2(a-6)2]- 


(275) 


Defining 


A  =  1  -  |p/|2|£|2, 

B  =  -2|p'|J(a •  £  -  |£|2), 

C  =  -|H2|a-6|2, 

X  =  A  +  Bx  +  Cx2, 


(276) 

(277) 

(278) 

(279) 


then  equation  (275)  can  be  written  , 


In  =  2/  dx  X- 


(280) 


or 
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I„  =  2 J  dx  (A  +  Bx  +  Cx2)-1. 


(281) 


Performing  the  integration,  then 


In 


2(B2  -  4AC)'1/2  in 


/2Cx  +  B  -  D  \  » 

UCx  +  B  +  D  /  o 


(282) 


where 


D  =  (B2  -  4AC),/2. 

Evaluating  the  upper  and  lower  limits,  equation  (282)  becomes 

In  =  2D'1  {in  [(2C  +  B  -  D)(2C  +  B  +  D)'1]  -  in  [(B  -  D)(B  +  D)'1]} 
or 


(283) 


(284) 


In  =  2D'1  in  [(2C  +  B  -  D)(B  +  D)(2C  +  B  +  D)'‘(B  -  D)'1]  (285) 


or 


In  =  2D'1  in  {[2C(B  +  D)  +  B2  -  D2][2C(B  -  D)  +  B2  -  D2]'1}. 
Substituting  equation  (283)  into  equation  (286), 

In  =  2D'1  in  {[2C(B  +  D)  +  4AC][2C(B  -  D)  +  4AC]'1} 
or 


(286) 


(287) 


In  =  2D'1  in  [(2A  +  B  +  D)(2A  +  B  -  D)-]. 
Using  equations  (276)  and  (277), 


(288) 
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2A  +  B  =  2(1  -  |p/|2|&|2)  -  2|p'|2(a-  6  -  |&|2) 


(289) 


or 


2A  +  B  =  2(1  -  |p'|2a •  6). 


(290) 


Substituting  equations  (244)  and  (245)  into  equation  (290),  then 


2A  +  B  =  2{l  -  |p'|2(2^  •  T)[ |C.|E/(|T|2  +  |^|2)]-‘}. 


(291) 


Clearly, 


2k  •  T  =  k  •  T  +  k  •  T. 


(292) 


By  equation  (243), 


■»  ♦  ♦ 
k  =  p  -  T. 


(293) 


Substituting  equation  (293)  into  the  first  term  on  the  right  side  of  equation  (292),  then 


21c  •  T  =  (p-T)-T  +  £  •  T,  (294) 

or 

2)c  •  T  =  p  -  T  -  |T|2  +  k  •  T.  (295) 


Substituting  equation  (243)  into  the  first  and  last  terms  on  the  right  of  equation  (294),  then 


2k-T  =  p-(p-k)-|T|2  +  Mp-  £) 


(296) 


or 


2^  •  T  =  |p|2  -  |t|2  -  |ic|2. 


(297) 
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Using  equations  (116)  and  (130), 


|p|2  =  (E'  +  co)2  -  m2.  (298) 

Substituting  equations  (21),  (236),  and  (298)  into  equation  (297),  then 

2k*  f  =  (E7  +  co)2  -  (E/2  -  |^|2)  -  |T I2  -  co2  (299) 

or 

2fc-T  =  |^|2  +  2E7co  -  |T|2.  (300) 

Then  substituting  equations  (21)  and  (300)  into  equation  (291), 

2A  +  B  =  2{l  -  frpflpf  +  2E'co  -  |t|2)[coE'(|t|2  +  fr|*)]-*}.  (301) 

Simplifying  equation  (301),  then 

2A  +  B  =  2(|t|2  -  |p'|2)(|T|2  +  Ipf^CL  +  |pf /coEO.  (302) 

Next,  to  evaluate  D  in  equation  (288),  one  first  evaluates  A,  B,  and  C.  Substituting  equation  (245)  in 
equation  (276), 

A  =  1  -  |f/|2/E/2  =  (E/2  -  |pf  )E'-2.  (303) 

Using  equation  (236),  equation  (303)  becomes 

A  =  m2/E/2.  (304) 

Further,  substituting  equations  (244)  and  (245)  into  equation  (277), 

B  =  -2|pf  {2 i  •  T[«E'(|f|2  +  |^|2)]-‘  -  E/_2}.  (305) 
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Also,  substituting  equations  (244)  and  (245)  into  equation  (278),  then 


C  =  -|p/|2[2T(|T|2  +  IH2 )-  -  £(a>E')"]2.  (306) 

Substituting  equations  (304)  to  (306)  into  equation  (283),  then 

D2  =  4|p/|4{(2£  •  t)2[oj2E/2(|t|2  +  |pf)*]" 

-  2  (2ic  •  T)[ojE/3(|t|2  +  |pf)]"}  +  4|p/|4/E/4 
+  4(m2|p/|2/E/2){4|f|2(|f|2  +  |^|2)-2 

-  2(2T  •  l<)[a)E/(|T|2  +  |^|2)h} 

+  4m2|p,|2/E/4.  (307) 

According  to  equation  (236), 

fpf  =  |p,|2(E/2  -  m2)  (308) 


and 


m2  =  E/2  -  fr\\  (309) 


Substituting  equations  (308)  and  (309)  into  the  second  and  third  terms,  respectively,  of  equation 
(307),  then 


D2  =  4|p'|4(2£-T)2  [co2E/2(|i|2  +  |^|2)2]-' 

+  4|*P'|2/E'2  +  4|p/|2{4|t|2(|t|2  +  |^|2)-2 
-  2(2T  •  ic)[coE/(|T|2  +  |p'|2)]-'} 

■  16|p/|4  |t|2[E/2(|t|2  +  l^l2)2]-1.  (310) 

Simplifying  equation  (310),  then 
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D2  =  4|p'|2(|t|2  +  |p'|2)'2[  |p/|2(2t  ■  T  )2(co2E/2)-‘ 

+  (|T|2  +  |H2)2/E'2-4|H2|T|2/E'2  +  4|T  |2 

-  2(2T  •  ic)(|T|2  +  |p/|2)(a)E/)"1]  (311) 

or 

D2  =  4|p'|2[a)2E/2(|f|2  +  |p/|2)2]~1[|p/|2(2£  •  T)2 
+  co2(|t|2  -  |p'|2)2  +  4a)2E/2|T|2 

-  2cuE/(2T  ■  &)(|t|2  +  |p/|2)] .  (312) 

Equivalently, 

D2  =  4|p'|2[o)2E/2(|f|2  +  |^|2)*]- 
{2i<  ■  t[  |p'|2(2ic  •  i)  -  2o)E/|t|2 
-  2a>E'|H2)  +  0)2(|T|2  -  |pf  )2 

+  4a)2E/2|f|2} .  (313) 

Next,  substituting  equation  (300)  into  equation  (313), 

D2  =  4|p'|2[o)2E/2(|f|2  +  |pf)2]-U  (314) 

where 


i d*  =  (Ip7!2  +  2E/o>  -  |T|2)(|p,|'t  +  2|J>fE/a>  -  |^|2|f|2 
-  2a>E'  |T |2  -  2a>E'  |^|2)  +  co2(|f|2  -  |^|2)2 
+  4a>2E/2|f|2. 


(315) 
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Multiplying  terms,  equation  (315)  becomes 


f„2  =  |£f  -  4coE/fPf  |t|2  -  2|pT|t|2  +  ZE'vfrr 

+  2coE'|T|4  +  |p'|2 |T|4  +  co2(|T|2  -  |^|2)2  (316) 

Dr 

=  rp/ii(ip/i4-2|p'i2ifi2  +  ifr  -  4coE^ifi2 

+  2E/co|p'|2)  +  2coE/|t|4  +  co2(|T|2  -  |pf)2  (317) 

or 

=  rP'|2(|T|2  -  |^|2)2  +  2coE/(|T|2  -  |^|2)2 

+  co2(|T|2  -  |pf)2  (318) 

or 

fD2  =  (|T|2-|pf)’(|pf  +  2coE'  +  CO2). 

Substituting  equation  (209)  into  equation  (319),  then 

f„2  =  (|t|2  -  |p'|2)2(E/2  -  m2  +  2coE/  +  co2) 

or 

fD2  =  (|T|2  -  frmw  +  co)2  -  m2].  (321) 


(319) 


(320) 


Next,  using  equation  (130)  in  equation  (321),  then 


Substituting  equation  (116)  into  equation  (322),  then 


h>2  =  (|T|2-|H2>2|p|2-  (323) 

Finally,  substituting  equation  (323)  into  equation  (314),  then 

D2  =  4|p|2|p'|2(|T|2  -  Ip*!2)2 

•  [cu2E/2(|f|2  +  fcf)2]-.  (324) 

Therefore, 

D  =  2|p| \j/\ (|T|2  -  |pf )  [o)E/(|f |2  +  |^|2)]-.  (325) 

Substituting  equations  (302)  and  (325)  into  equation  (288),  then 

In  =  £*>E/(|T|2  +  |H2)[|p|  |H(|i|2  -  |p|2)]-*  in  A,„  (326) 

where 


An  =  {2(|T|2  -  |pf)(|T|2  +  |P'|2)-,(1  +  |p'|2/wE/) 

-  2|p|  |p'|(|T|2  -  |p'|2)[a>E/(|f|2  +  |H2)]-*} 

•  {2(|T|2  -  |p'|2)(|t|2  +  |pf)-(l  +  |p'|2/coE/) 
+  2|p|  |p'|(|f |2  -  |p/|2)[coE/(|t|2  +  fpf)]-1}" 


or 


A„  =  (1  +  If/IVcoE2  -  |p|  |p'|/coE/) 

•  (1  +  |p,|2/cpE/  +  |p|  |p'|/coE/)"1 


(327) 


(328) 


or 
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A,,  =  (coF  +  IH2-  IpUHXqiE'  +  +  IpIIp'I )-*.  (329) 

Substituting  for  co  in  equation  (329)  with  equation  (130),  then 

A„  =  [(E-E')E'  +  fr\'-\m] 

•  [(E-E')E'  +  fr\*  +  IpIIH]-  (330) 

or 

An  =  (EE'-F2  +  |H2- IpIIp'I) 

•  (EE'-E/2  +  Ip'I2  +  IpIIp'I)'1.  (331) 

Using  equation  (209)  in  equation  (331),  then 

An  =  (EF  -  m2  -  |p|  |^|)(EF  -  m2  +  |p|  |^|)-‘.  (332) 

Finally,  substituting  equation  (332)  into  equation  (326),  then 

In  =  coF(|t|2  +  fpf  )[|pI  |p'|(|T|2  -  fp'|2)]-‘L,  (333) 

where 

L  =  in  [(EF  -  m2  +  |p|  |^|)(EF  -  m2  -  |p|  |p'|)-‘].  (334) 

Proceeding  with  the  evaluation  of  I12  by  using  equation  (252),  then 

I12  =  f  (d2Q//2n)(l  -  ^  •  i)-‘(l  -  P'  *  i)-2.  (335) 

Using  the  identity  equation  (266),  one  first  notes  that 
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d/dp{ap)-1  =  a-ld/dp(fi-1)  =  -a"1/?"2. 


(336) 


and  therefore 


or 


a~lp~2  =  -d/dp{ap)-'  =  -d/dp  f  dx[ax  +  jj(  1  -  x)]'2 

0 


(337) 


a-ip-i  =  j-  dx  2(1  -  x)[ax  +  P(1  -  x)]-3. 
0 


(338) 


Using  the  identity  equation  (338)  in  the  integrand  of  equation  (335),  then 


1.2  =  /'  dx/(d2QV2n)2(l  -  x)[(l  -  p>  •  a)x  +  (1  -  p'  •  £)(1  -  x)]'3  (339) 

0 


or 


1.2  =  /  dx  2(1  -  x)/(d2QV2n)[l  -  p>  -g]'3, 
0 


(340) 


where  g  is  given  by  equation  (271). 


Defining, 


I.2»(x)  =  J(d2QV2n)[l  -  ^  •  g]'3, 


(341) 


then  equation  (340)  is  written  as 


Il2  =  /  dx  2(1  -  x)Ii2o(x). 
0 


(342) 


Again  choosing  the  polar  axis  of  integration  to  lie  along  the  vector  g,  then  equation  (341)  becomes 


I.2o(x)  =  /"  d6'  sin  61  [l-  J p7 1  jgjcos  0']-3. 
0 


(343) 
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Changing  variables  of  integration  according  to 


y  =  Ip'IIsI  cos  9', 


(344) 


then  equation  (343)  becomes 


Ipllgl 

If2o(x)  =  (fp'Hil)-  /  _  .  dy(l-y)- 
-Ipllgl 


(345) 


or 


Ipllgl 

Ino(x)  =  (2|f/|  |gj)-1(l  -  y)~2  | 


-Ipllgl 


(346) 


or 


I120(x)  =  2(1  HM2)-2. 


(347) 


Substituting  equation  (347)  into  equation  (342),  then 

I12  =  4  /  dx  (1  -  x)X"2,  (348) 

0 

where  X  is  defined  by  equations  (276)  to  (279).  The  reduction  to  this  form  is  identical  to  that  of  equa¬ 
tion  (273)  to  equation  (281).  Equation  (348)  is  rewritten  as 


I12  —  4li2i  -  4li22,  (349) 


where 


(350) 


1 122  —  /  dx  xX"2. 
0 


(351) 
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Performing  the  integral,  equation  (350)  becomes 


Im  =  (2Cx  +  B)(4AC  -  B2)-(Cx2  +  Bx  +  A)-‘  I 

0 

"b  2C(4AC  -  B2)  1  /  dx  (Cx2  +  Bx  +  A)"1 


(352) 


or  evaluating  the  first  term  and  using  equations  (281)  and  (283),  then 


1.2.  =  -(2C  +  B)[D2(C  +'B  +  A)]-'  +  B(D2A)-'  -  CD'2!,,. 


(353) 


Combining  terms,  equation  (353)  becomes 


Im  =  (B2  -  2AC  +  BC)[AD2(A  +  B  +  C)]'1  -  CD'2!,,, 


(354) 


or  using  equation  (283),  then 


1.2.  =  (D2  +  2AC  +  BC)[AD2(A  +  B  +  C)]-  -  CD'2!,,. 


(355) 


Next,  evaluating  equation  (351), 


I122  — (2A  +  Bx)(4AC  -  B2)  ‘(A  +  Bx  +  Cx2)'1 
-  B(4AC  -  B2)'1  /*  dx  (A  +  Bx  +  Cx2)'1. 


(356) 


Again  using  equations  (281)  and  (283),  equation  (356)  becomes 


1.22  =  -(B  +  2C)D'2(A  +  B  +  C)'1  +  (B/2D2)I„. 


(357) 


Next,  substituting  equations  (355)  and  (357)  into  equation  (349),  then 


I,2  =  4(D2  +  2AC  +  BC)[AD2(A  +  B  +  C)]'1  -  4CD'2I„ 

+  4(B  +  2C)D'2(A  +  B  +  C)'1  -  2BD'2I„.  (358) 
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Combining  terms,  then  equation  (358)  becomes 


I12  =  4(D2  +  4AC  +  BC  +  AB)[AD2(A  +  B  +  C)]-  -  2  (B  +  2C)D-2I„. 


(359) 


Using  equation  (283)  in  the  numerator  of  the  first  term,  then 


I12  =  4(B2  +  BC  +  AB)[AD2(A  +  B  +  C)]-  -  2(B  +  2C)D-2I,, 


(360) 


I12  =  4B(AD2)-  -  2(B  +  2C)D-2IU 


(361) 


Next, 


I2,  =  /  (d2£2//2n)(l  -  p'  •  a)‘2(l  -  p'  •  fi)-. 


(362) 


Using  the  identity  equation  (266),  one  notes  that 


a ~2p~'  =  -d/d a{afiY'  =  -d  /  d a  J  dx  [ax  +  /3(1  -  x)]"2 


(363) 


or 


a'2/ 3-'  =  2  /  dx  x[ax  +  ft(l  -  x)]'\ 


(364) 


Rewriting  the  integrand  of  equation  (362)  with  the  identity  equation  (364)  and  the  definition  equation 
(271),  then 


I21  =  2  /  dx  x/  (d2S2'/2n)(l  -  p'  •  ^)- 


(365) 


Substituting  equation  (341)  into  equation  (365),  then 
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(366) 


111  —  2  f  dx  xIi2o(x). 
0 

Then  substituting  equation  (347)  into  equation  (366), 


l2.=4/  dx  x(l  -  Ip'Hgl2)'2. 

0 

Again  recalling  the  reduction  of  equation  (273)  to  equation  (281)  and  using  equations 
then  equation  (367)  becomes 

I21  =  4  /  dx  xX'2. 


Comparing  equations  (368)  and  (351),  and 


III  —  41]  22  • 


Finally,  substituting  equation  (357)  into  equation  (369),  then 


In  =  -4(B  +  2C)D'2(A  +  B  +  C)"  +  2BD'2IU. 


Using  equation  (252),  I22  is  given  by 


In  =  J  (d2S2'/2n)(l  -  }/  •  a)-2(l  -  ■  i)'\ 


One  notes  using  equation  (266)  that 

a  2P  2  =  (d  /  da)(d  / d p)(ap)~'  =  d1/  d a  d p  /  dx[ox  +  p(  1  -  x)]-2 


or 


a  2p  1 - 23  /  da  /  dx  (1  -  x)[ax  +  p(  1  -  x)]  3, 


or 


(367) 

(276)  to  (279), 

(368) 

(369) 

(370) 

(371) 

(372) 

(373) 
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(374) 


a  2ft  2  =  6  f  dx  x(l  -  x)[crx  -f  p(l  -  x)]  4. 
o 


Using  the  identity  equation  (374)  in  equation  (371),  then 


I22  =  6  /  dx  x(l  -  x)I22i(x), 
o 


(375) 


where,  analogously. 


I221  (x)  =  J  (d2Q'/2n)(l  -  p7  •  g)~4. 

Again  choosing  the  polar  axis  of  integration  along  the  vector  g,  then 

I221U)  =  /  dQ'  sin  O'  (1  -  |p'|  |g|  cos  Q'Y 4, 
0 

and  making  the  substitution  with  equation  (344),  then 

fp'llgl 

I22 1  (x)  =  |p'|_I  IsT1  f  dy(l-y)-4. 
-Ip'llsl 

Performing  the  integration,  then  equation  (378)  becomes 


(376) 


(377) 


(378) 


Wx)  =  (3|p'||g|)-[(l-|p'l|g|)-3 

-(1  +  IHIM. 


(379) 


Simplifying  equation  (379),  then 


Wx)  =  (2/3)(3  +  |^|2|i|I)(l-|p'l2|g|2)'3- 


(380) 


Then  substituting  equation  (380)  into  equation  (375), 


I22  =  6  /  dx  x(l  -  x)(2/3)[3  +  (1  -  X)]X'3, 
0 


(381) 
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where  equation  (279)  has  again  been  used.  Rewriting  equation  (381),  then 


I22  —  I6I221  -  I6I222  ~  4I122  "b  4l223/  (382) 

where  equation  (351)  has  been  used,  and 

1221  =  /  dx  xX“3,  (383) 

o 

1222  =  /  dx  x2X  3,  (384) 

0 

1223  =  /  dx  x2X~2.  (385) 

0 

Evaluating  the  integral  of  equation  (383), 

I221  =  -(2A  +  Bx)[2(4AC  -  B2)(A  +  Bx  +  Cx2)2]'1  |  *  -  3B[2(4AC  -  B2)]-1  /  dx  X'2.  (386) 

0  o 

Evaluating  the  first  term,  using  equation  (283),  and  substituting  equation  (350),  then 

I221  =  (2A  +  B)[2D2(A  +  B  +  C)2]’1  -  D-2A-  +  (3/2)BD-2I12l.  (387) 

Substituting  equation  (355)  into  equation  (387),  then 

I221  =  (2A  +  B)[2D2(A  +  B  +  C)2]'1  -  D'2A-'  +  (3/2)BD-2 

•  [(D2  +  2AC  +  BC)A~‘D~2(A  +  B  +  C)'1  -  CD^hJ.  (388) 

Combining  terms,  then  equation  (388)  becomes 

I22.  =  [(2A  +  B)D2A  -  2D2(A  +  B  +  C)2  +  3B(D2  +  2AC  +  BC)(A  +  B  +  C)] 

•  [2D4A(A  +  B  +  C)2]-  -  (3/2)BCD-4In.  (389) 

Next,  evaluating  I222  of  equation  (384),  then 
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f 


l222  =  -x[3C(A  +  Bx  +  Cx2)2]'1  |  -  (B/3C)  /  dx  x(A  +  Bx  +  Cx2)'3 

o  o 

+  (A/3C)  /  dx  [A  +  Bx  +  Cx2]-3. 

0 


(390) 


Evaluating  the  first  term  of  equation  (390)  and  using  equation  (383),  then 


I222  =  -[3C(A  +  B  +  C)2]-1  -  (B/3QI221  +  (A/3C)I222i,  (391) 


where 


I2221  =  /  dx  (A  +  Bx  +  Cx2)-3. 
0 


(392) 


Evaluating  equation  (392),  then 


I2221  =  (2Cx  +  B)[2(4AC  -  B2)(A  +  Bx  +  Cx2)2]'1  |  ^  +  6C[2(4AC  -  B2)]'1  Jf'  dx  X'2.  (393) 


Evaluating  the  first  term  and  using  equations  (350)  and  (355),  then  equation  (393)  becomes 


I2221  =  (2C  +  B)[2(4AC  -  B2)(A  +  B  +  C)2]-  -  B[2(4AC  -  B2)A2]- 

-  3CD'2{(D2  +  2AC  +  BC)[AD2(A  +  B  +  C)]"1  -  CD-2I„}.  (394) 

Simplifying  equation  (394)  and  using  equation  (283),  then 

T  =  2C  +  B  ,  B  3C(D2  +  2AC  +  BC)  ,  3C*_  ,  f,p(-x 

12221  '  2D2(A  +  B  +  C)2  2D2A2  "  AD4(A  +  B  +  C)  D4 


Then  by  substituting  equations  (389)  and  (395)  into  equation  (391), 

1  B  ("  (2A  +  B)D2A-2D2(A  +  B  +  C)2  +  3B(D2  +  2AC  +  BC)(A  +  B  +  C) 

1222  ~3C(A  +  B  +  C)2  "  3C  I  2D4A(A  +  B  +  C)2 

_3  BC  T  I  +  A  f  2C  +  B  ,  B  3C(D2  +  2AC  +  BC)  3C2  ,  1 
"2  D4  "J  3C  L’  2D2(A  +  B  +  C)2  2D2 A2  *  AD4(A  +  B  +  C)  D4  "J ’ 

(396) 


Next,  by  evaluating  I223  of  equation  (385),  then 


64 


X 


(397) 


I223  —  ■ 


C(A  +  Bx  +  Cx2) 


+ 


A 

C 


/  dx  X- 


Evaluating  the  first  term  and  using  equations  (350)  and  (355),  then  equation  (397)  becomes 


I 


223 


1 _ 

C(A  +  B  +  C) 


D2  +  2AC  +  BC  C_  1 
AD2(A  +  B  +  C)  '  D2  'j 


Simplifying  equation  (398)  and  using  equation  (283),  then 


(398) 


2A  +  B  _  A. 
D2(A  +  B  +  C)  D2 


(399) 


Finally,  substituting  equations  (389),  (396),  (357),  and  (399)  into  equation  (382),  then 

I22  =  16(1  +  B/3C){[2D4A(A  +  B  +  C)2]-'[(2A  +  B)D2A  -  2D2(A  +  B  +  C)2 

+  3B(D2  +  2AC  +  BC)(A  +  B  +  C)]  -  (3BC/2D4)I„}  -  16(-[3C(A  +  B  +  C)2]'1 
+  (A/3C){-(2C  +  B)[2D2(A  +  B  +  C)2]'1  +  B/2D2A2  -  3C(D2  +  2AC  +  BC) 

•  [AD4(A  +  B  +  C)]-1  +  (3C2/D4)I„})  -  4{-(B  +  2C)[D2(A  +  B  +  C)]’1 
+  (B/2D2)Ii,}  +  4{(2A  +  B)[D2(A  +  B  +  C)]’1  -  (A/D2)!,,}.  (400) 


Combining  terms,  equation  (400)  becomes 


I22  =  [i6(l  +  B/3C)(-3BC/2D4)  -  (16A/3C)(3C2/D4)  -  2B/D2  -  4A/D2]l„  +  16(1  +  B/3C) 

•  [2D4A(A  +  B  +  C)2]-‘[(2A  +  B)D2A  -  2D2(A  +  B  +  C)2  +  3B(D2  +  2AC  +  BC) 

•  (A  +  B  +  C)]  -  16(-[3C(A  +  B  +  C)2]-'  +  (A/3C){-(2C  +  B)[2D2(A  +  B  +  C)2]- 

+  B/2D2A2  -  3C(D2  +  2AC  +  BC)[AD4(A  +  B  +  C)]'1})  +  8/D2.  (401) 


Next,  using  equation  (252),  I2-i  is  given  by 


I2-,  =  /  (d2Q72ir)  (1  -  p'  ■  a)-2(l  -!/■{,). 


(402) 
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Define  a  unit  vector  n  perpendicular  to  a  and  6,  namely. 


n  -  a  X  b/  lallb  I 


(403) 


and  choose  a  vector  nj^  perpendicular  to  a  and  n,  namely. 


n_|_  =  a  X  n  /\a\. 


(404) 


The  vectors  n,  n±,  and  a  define  orthogonal  coordinate  axes  for  the  integration.  The  polar  axis  is  taken 
along  the  vector  a.  The  vector  p'  can  be  written  in  terms  of  its  projections  along  n,  n^,  and  a  as 
follows: 


P'  =  (p7  *  a/|a|2)a  +  p'  •  n  n 
+  P'  '  n±  . 


(405) 


Then 


p-6  =  p-aa-6/|ap 

+  p/  *  n  b  *  n  +  p—  n_j_  b  *  nj_ . 


(406) 


Using  equation  (403),  one  notes  that 


b  •  n  =  0. 


(407) 


Choosing  the  polar  axis  of  integration  to  lie  along  the  vector  a  and  the  azimuthal  angle  to  be  with 
respect  to  n±/  then 


p1 '  n_|_  =  |p'|  sin  9'  cos  . 


(408) 


Substituting  equations  (407)  and  (408)  into  equation  (406),  then 
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P'*  b  =  (P'  *  a/|a|2)a  *  b 

+  b  •  nj_  Ip7!  sin  9*  cos  <j>; . 


(409) 


Using  equation  (409)  in  equation  (402),  then 


I2-i  =  /  (2tt)  1  sin  9'  dQ/  d<j>'(l  -  p7  •  a)  2 

*  (1  -  p/  *  a  a  *  b  | a | ~2  -  b  *  nj_  |p;|  sin  9/  COS  <j>;) 

=  J(2n)~l  sin  9'  dO'  d+'(l  -  p'  •  a)‘2  (1  -  ^  •  a  a  ■  b|a|'2)  +  I2_n,  (410) 

where 

I2-11  =  -J(2n)-1  sin  d0;  d^l  -  p'  •  a)-2  b  *  njjp/l  sin  cos  <j>;.  (411) 


But  clearly 


I2-11  =  -  J(2n)  1  sin2  9f  (1  -  |p/|  lajcosGO^dG'  i  •  Aj_  Ip7!  /  cos  <t>' d<j>' 

0 


(412) 


or 


I2-11  —  0. 


(413) 


Substituting  equation  (413)  into  equation  (410),  then 

I2-.  =  /( 2n)-  d2Q'  (1  -  p'  •  a)'2(l  -  y  •  a  a  ■  6|a|-2).  (414) 


Rewriting  equation  (414)  and  then  using  the  above  coordinate  frame, 


I2-,  =  J (2 n)“  d2Q'  (1  -  ■  a)’2  -  a  •  6|a|-2 

•  J (2 n)"1  d2Q2  p/  •  a(l  -  p'  ‘  a)-2 


(415) 
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or 


[2-,  =  / (27 z)-‘  d2Q'  (1  -  •  a)'2  -  a  •  b|a|“2 

•  /( 277)*'  d2Q7  (1  -  p'  ■  a)"2[l  -  (1  -  p'  •  a)] 


(416) 


or 


I2-.  =  (1  -  a  ■  b|a|-2)  / (27l)-*  d2Q'  (1  -  \t  •  a)’2 
+  a  ■  b|a|‘2  / (27i)-‘  d2£2/  (1  -  p'  •  a)"‘. 

Comparing  equation  (417)  with  equations  (263)  and  (261),  then 


I2-,  =  (1  -  a  •  6|a|-2)I20  +  a  •  £|a|~2I10. 


(418) 


Finally,  substituting  equations  (262)  and  (265)  into  equation  (418),  then 


I2-,  —  a  •  6|a|~2(|p/|  |a|)-’  in  [(1  +  |^|  |a|) 

’(I-  |p'l|a|)“]  +  2(1  -  a  •  b|a|~2)(l  -  |p/|2|a|2)-1. 


(419) 


To  proceed  then,  comparing  equations  (238)  and  (239),  then  equation  (238)  may  be  rewritten  as 


+  2(E2  +  E/2)]T „  -  x-'To,  -  [2x-‘(w  -  |p|  cos  9)+  |p|2x'2  sin2  9  -  l]l10 
-  4m2E2T22  +  xT^E'E2*  +  2w2|p|2  sin2  9  +  2cu2|^|2  -  2cu2E/2  -  2(E2  +  E/2)(|p|2 
+  | P'12  +  to2  -  2|p|w  cos  9  +  2coE'  -  2|p|E'  cos  0)]l21  +  [4x'‘(E2  +  E/2) 

•  (w  -  |p|  cos  9)  -  4E2  +  4x"2|p|2E/2  sin2  9  +  4E/coV‘]T20  -  2o>VT2.,. 

Substituting  equation  (251)  into  equation  (420),  then 
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I  =  m2(|T|2  +  |p'|2)'1E/"2Ii2  +  [-2E'x  +  |p|2  +  fcf  +  cp2 


-  2|p|co  cos  0  +  2cpE'  -  2 |p|  cos  0  V  +  2(E2  +  E/2)](|t|2  +  |^|2)-(E/x)-,Iu 

-  jc-E'-Io,  -  [2x“‘(tp  -  |p|  cos  0)  +  |p|2x"2sin2  0  -  l](|f|2  +  |p,|2)"1Iio 

-  4m2E2(|T|2  +  |^|2)-2E/'2I22  +  jf^E'E2*  +  2cp2|p|2  sin2  0  +  2cp2|[>'|2 

-  2cp2E/2  -  2(E2  +  E/2)(|p|2  +  |^|2  +  cp2  -  2|p|cp  cos  0  +  2wE'  -  2|p|E2  cos  0)] 

•  (|T|2  +  |^|2)-E-‘I21  +  [4x-(E2  +  E/2)(cp  -  |p|  cos  0)  -  4E2 

+  4|p|2E/2x'2  sin2  0  +  4E/cp2x~1](|t|2  +  |p/|2)-2I20  -  2cp2x"1(|t|2 
+  |p/|2)_2E/I2_i.  (421) 

Substituting  equations  (361),  (260),  (262),  (401),  (370),  (265),  and  (419)  into  equation  (421),  then 

I  =  m2(|T|2  +  |pf  )-‘E/'2[4B/AD2  -  2(B  +  2C)D-2In]  +  (|t|2 
+  |f/|2)'lE/“‘x'1[-2E/x  +  |p|2  +  |^|2  +  cp2  -  2|p|cp  cos  0 
+  2cpE2  -  2|p|E2  cos  0  +  2E2  +  2E/2]l„  -  x-,E/-,|p'|-,|6|-1 

•  in  [(1  +  |p'||£|)(l  -  fp'Hbl)-]  -  [2x-(cp  -  |p|  cos  0) 

+  |p|V2  sin2  0  -  l](|T|2  +  |i>f)-^|-‘fa|-> 

•  in  [(1  +  |p'||a|)(l  -  |^|  |a|)-]  -  4m2E2E/'2(|f|2  +  |^|2)-2 

•  [16(1  +  B/3C)(-3BC/2D4)  -  16ACD-4  -  2BD~2  -  4AD‘2]lu  -  4m2E2E/'2(|T|2  +  fpf)-2 

•  (16(1  +  B/3C)[(2A  +  B)D2A  -  2D2(A  +  B  +  C)2 

+  3B(D2  +  2AC  +  BC)(A  +  B  +  C)][2D4A(A  +  B  +  C)2]’1 
+  16  [3C(A  +  B  +  C)2]-1  -  (16A/3C){(2C  +  B)[-2D2(A  +  B  +  C)2]’1 
+  B/2D2A2  -  3C(D2  +  2AC  +  BC)[AD4(A  +  B  +  C)]'1}  +  8D-2) 

+  xT^E'E2*  +  2cp2|p|2  sin2  0  +  2cp2|^|2  -  2cp2E/2  -  2(E2  +  E/2)(|p|2 
+  Ip2!2  +  cp2  -  2|p|cp  cos  0  +  2cpE'  -  2|p|E/  cos  0)] 

■  (|t|2  +  rP/|2)-2E/“{-4(B  +  2C)[D2(A  +  B  +  C)]1 
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(422) 


+  2BD  2I,,}  +  [4(E2  +  E/2)x  ‘(cp  -  |p|  cos  9)  -  4E2  +  4|p|2E/2x'2  sin2  6 
+  4E/cp2x'1](|t|2  +  |p/|2)_22(l  -  |pf  |a|2)- 
-  2<p2x'*(|T|2  +  |p,|2)"2E/{a  •  6|a|-3|J/|- 
•  in  [(1  +  |p/|  |a|)(l  -  |p/|  |a|)_1  ] 

+  2(1  -  a  •  £|a|"2)(l  -  |p/|2|a|2)-'}  . 

Separating  terms  containing  IM,  then  equation  (422)  may  be  rewritten  in  the  form 

I  =  I.  +  I»,  (423) 

where 

l  =  {-2(B  +  2C)m2(|f  |2  +  |J>f)-*E'-2D-2  +  (|f|2  +  fpf  J-x-E'-^E'x 

+  IpI2  +  Ip2!2  +  co2  -  2|p|co  cos  6  +  2cpE'  -  2|p|  cos  9  E'  +  2E2  +  2E/2) 

-  4m2E2E/‘2(|f  |2  +  |p'|2)_2[-24(l  +  B/3C)BCD-4  -  16ACD'4  -  2BD'2  -  4 AD'2] 

+  2Bx"1D"2E/~'(|T|2  +  l^^'dSE'E^  +  2cp2|p|2  sin2  9  +  2cp2|p/|2  -  2cp2E/2  -  2(E2 
+  E/2)(|p|2  +  |pf  +  cp2  -  2|p|cp  cos  9  +  2cpE2  -  2|p|E'cos  0)]}ln,  (424) 

It  =  m2(|T|2  +  |p,|2)'1E/'2(4B/AD2)  -  (xFIp'IISI)'* 

,1 

•in  [(1  +  |p'||b|)(l  -  |p'||fi|)'1]  -  [2x'*(cp  -  |p|  cos  9) 

+  |p|2x'2  sin2  9  -  1](|T|2  +  |^|2)-1|p2|-|a|- 

•in  [(1  +  |p/| |a|)(l  -  |p'| |a| )_1  ]  -  4m2E2E/'2(|f|2 

+  |p/|2)'2{l6[l  +  (B/3C)](2D4A)-'(A  +  B  +  C)'2[(2A  +  B)D2A 

-  2D2(A  +  B  +  C)2  +  3B(D2  +  2AC  +  BC)(A  +  B  +  C)] 

+  (16/3C)(A  +  B  +  C)-2  -  (16A/3C)[(2C  +  B)(-2D2)'*(A  +  B  +  C)'2 
+  (B/2D2A2)  -  3CA''D'4(A  +  B  +  C)''(D2  +  2AC  +  BC)]  +  8D'2} 

+  x-D-dTl2  +  |pf)'2[-4(B  +  2C)D'2(A  +  B  +  C)-][8E'E2x  +  2cu2|p|2  sin2  9 
+  2oj2|p/|2  -  2<p2E/2  -  2(E2  +  E/2)(|p|2  +  |p'|2  +  cp2  -  2|p|cp  cos  9 
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+  2ajE7  -  2|p|E'  cos  0)]  +  [4x  ’(E2  +  E/2)(a)  -  |p|  cos  0)  -  4E2  +  4x‘2|p|2E'2  sin2  0 
+  4E/x~'a>2](|T|2  +  |p'|2)*22(l  -  |p/|2|a|2)”1 
-  2a>2x"'(|T|2  +  |p/|2)’2  E'{a  •  blp'l  1  |a|  3 

•  in  [(1  +  IH  |a|)(l  -  IVI  |a|)-‘]  +  2(1  -  a  •  6|a|-2) 

•  (1-  fpffal2)-1}-  (425) 

Substituting  equation  (333)  into  equation  (424),  then 

I«  =  {Li  +  x-'W-'irlWx  +  |p|2  +  fpf  +  co2  -  2|p|co  cos  6  +  2coE'  -  2|p|E'  cos  6 
+  2E2  +  2E/2)  -  4m2E2E/_2(|f|2  +  |pf  )"lFal  +  2Bx-,D-2E/“1(|t|2  +  Ipf)"1 
•  [8E/E2x  +  2co2|p|2  sin2  6  +  2cd2\p'\2  -  2cd2E'2  -  2(E2  +  E/2)(|p|2  +  |pf  +  a;2 
-  2|p|co  cos  6  +  2wE'  -  2|p|E7  cos  6)]} E'cofpt1  |p|_1(T2  -  IpfPL,  (426) 

where 

Ial  =  -2(B  +  2C)m2E/_2D~2,  (427) 

Fo1  =  -24(1  +  B/3C)BCD'4  -  16ACD"4  -  2BD“2  -  4 AD-2.  (428) 

To  reduce  equation  (427),  one  first  substitutes  equation  (243)  into  equation  (305),  and  using  equa¬ 
tions  (191)  and  (41),  then 

B  =  -2|p/|2[(2|p|gl>  cos  6  -  2<jo2)co“1E/"1(|t|2  +  |p' |2)-1  -  E/_2)].  (429) 

Equivalently, 

B  =  -2|p/|2(2E/ |p|co  cos  0  -  2E/co2  -  co|T|2  -  co|p'|2)co"1E/“2 

*  (|t|2  +  |pfr.  (430) 


71 


Substituting  equation  (243)  into  the  numerator  of  equation  (430)  and  using  equation  (41),  then 


B - 2|p,|2cj"‘E/  2(|f|2  +  Ip7!2)  ‘(ZE'IpIo)  cos  9  -  2E'w2  -  co|p|2  -  co3 

+  2|p|w2  cos  6  -  o)|p'|2).  (431) 

Substituting  equation  (130)  into  equation  (431),  then  equation  (429)  becomes 

B  =  -2|p/|2E/_2(|T|z  +  |f/|2)-‘(2|p|E  cos  6  -  2Eco  -  |p|2  +  co2  -  |p7|2).  (432) 

Also,  by  substituting  equation  (130)  into  equation  (209)  and  using  equation  (116),  then 

|pf  =  (E  -  a))2  -  m2 

=  E2  -  2Egd  4*  co2  -  m2 

=  |p|2  -  2Eco  +  w2.  (433) 

Substituting  equation  (433)  into  equation  (432),  then 

B  =  -2|p'|2E/"2(|T|2  +  |p/|2)'1(2|p|E  cos  9  -  2|p|2).  (434) 

Furthermore,  using  equation  (43)  in  equation  (434),  then 

B  =  -2|p/|2E/-2(|t|2  +  |p/|2)'1[2E(E  -  x)  -  2|p|2],  (435) 

or  again  using  equation  (116),  then 

B  =  -4|p/|2E/'2(|f|2  +  |p'|2)-‘(m2  -  xE).  (436) 

Further,  using  equation  (41)  in  equation  (306),  then 

C  =  -|p,I2[4|t|2  (|T|2  +  |J/|2)-2  -2(2t-T) 

•  w-'E'-'dTl1  +  Ip7!2)’1  +  E2'2].  (437) 
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Equation  (437)  may  be  rewritten 


C  =  -|p'|Iw“E/'2(|f|2  +  |pf  H4wE'2  |t|2  -  2(21c  •  T) 

•  E'(|T|2  +  fpf )  +  w(|t|2  +  |pf  )2]  (438) 

or 

C  =  -|p'|2E/*2(|f|I  +  |p'|2)'2FC/  (439) 

where 

F c  =  4E/2|f|2  -  2{2l  •  T)aT’E'(|T|2  +  fp'l2) 

+  (|T|2  +  fp'l2)2.  (440) 

Equation  (440)  can  be  written 

Fc  =  (2t  •  ToT'E'  -  |T|2  -  |p'|2)2  -  4(t  •  T)2E/2oT2 

+  4E/2|T|2.  (441) 

Substituting  equation  (243)  and  using  equation  (191)  in  equation  (441),  then 

Fc  =  [2aT‘(u>|p|  cos  6  -  co2)E/  -  (|p|2  -  2|p|co  cos  6  +  u>2)  -  |p'|2]2 

-  4E/2[w'2(a)|p|  cos  6  -  to2)2  -  |p|2  +  2|p|«  cos  6  -  co2].  (442) 

Substituting  equation  (130)  into  equation  (442),  then 

Fc  =  [2E(|p|  cos  0  -  to)  +  co2  -  |p|2  -  |p'|2]2  +  4E/2|p|2(l  -  cos2  0).  (443) 

Substituting  equations  (116)  and  (209)  and  using  equations  (130)  and  (213)  in  equation  (443),  then 

Fc  =  (2E|p|  cos  0  -  2E2  +  2m2)2  +  4E/2|p|2  sin2  0.  (444) 
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Substituting  equation  (43)  into  equation  (444),  then 


Fc  —  4[(m2  -  jcE)2  +  E/2|p[2  sin2  0],  (445) 

Then,  substituting  equation  (445)  into  equation  (439), 

C  =  -4|pfE'-2(|T|2  +  |p'|2)~2[(m2  -  xE)2  +  E/2|p|2  sin2  0],  (446) 

Also,  using  equations  (243),  (191),  and  (41), 

|f|2  -  Ip7!2  =  |p|2  -2|p|co  cos  0  +  co2  -  |p,|2.  (447) 

Substituting  equations  (116)  and  (209)  into  equation  (447),  substituting  equation  (130),  and  using 
equation  (43),  then 


|T|2  -  |p'|2  =  2 cox.  (448) 

Substituting  equation  (448)  into  equation  (325),  then 

D  =  4|p|  |p'|kE'-,(|T|2  +  |^|2)-*.  (449) 

Next,  by  substituting  equations  (436),  (446),  and  (449)  into  equation  (427),  then 

Io1  =  m2(2E'2|p|2x2)-[(m2  -  xE)(|T|2  +  |pf)  +  2(m2  -  xE)2  +  2E/2  |p|2  sin2  0],  (450) 

Using  equations  (243),  (191),  and  (43) 

Jt|2  +  Ip7!2  =  |p|2  -  2co(E  -  x)  +  co2  +  |f/|2.  (451) 

Also,  using  equations  (213)  and  (43), 

|p|2  sin2  0  =  |p|2  -  E2  +  2Ex  -  x2.  (452) 
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Then,  substituting  equations  (451)  and  (452)  into  equation  (450), 


I«i  =  m2(2E/2|p|2x2)  *{(m2  -  xE)[|p|2  -  2co(E  -  x)  +  co2  +  IpI2] 
+  2(m2  -  xE)2  -  2E/2|p|2  -  2E/2E2  +  4E/2Ex  -  2E/2x2}. 


(453) 


Substituting  equations  (116)  and  (209)  and  using  equation  (130)  in  equation  (453),  then 


Lt  =  m2(2E/2|p|2x2)-‘[(m2  -  xE)(2E'2  -  2E'x)  -  2E'2m2  +  4E'2Ex  -  2E'V]  (454) 


or 


I„,  =  m2(2E/2|p|2x2)”‘(-2m2E/x  -  2xEE'2  +  2x2EE'  +  4E'2Ex  -  2E/2x2). 


(455) 


Using  equations  (116)  and  (130), 


-2m2E7x  =  -2E2E7x  4-  2|p|2E7x 

=  -2E(E7  4-  co)E7x  +  2|p|2E7x. 


(456) 


Also,  using  equation  (130), 


2x2EE7  =  2x2(E7  4-  co)E7. 


(457) 


Substituting  equations  (456)  and  (457)  into  equation  (455),  then 


I«i  =  m2(2E/2|p|2x2)  ^^EcoE'x  4-  2|p|2E7x  4-  2x2E7co), 


(458) 


or 


Ial  =  m2E7  *|p|  2x~1(-Eco  4-  |p|2  4-  xco) 


(459) 


or 
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-  m2a)EE,',|p|-V  +  m2wE/‘,|p|-2. 

Equation  (428)  can  be  written  as 

F„i  =  -2D-2[4D-2(3BC  +  B2  +  2AC)  +  B  +  2A], 

Using  equation  (283)  in  equation  (461),  then 

F„,  =  -2D‘2[4  +  12Fo11  +  B  +  2A], 

where 

F„n  =  (2A  +  B)CD-2. 

Substituting  equations  (304),  (449),  (436),  and  (446)  into  equation  (463),  then 

F°n  =  -[(m2  -  xE)2  +  E/2|p|2  sin2  0]  [2|p|2x2E/2(|f|2  +  |pf  )]-*F0ll  , 

where 

Fain  =  m2|f|2  +  m2|^|2-2m2|^|2  +  2xE|^|2. 

Substituting  equation  (451)  into  equation  (465),  then 

F«iii  =  m2|p|2  -  2ojm2(E  -  x)  +  m2ca2  -  m2|p/|2  +  2xE|p/|2. 
Substituting  equations  (116)  and  (209)  into  equation  (466),  then 

F.iii  =  m2E2  -  m2E/2  +  m2co2  -  2Ewm2  +  2x(E|p'|2  +  cam2). 
Substituting  equation  (130)  into  equation  (467), 

F„in  =  2x(E|p'|2  +  cj2m2). 


(460) 


(461) 


(462) 


(463) 


(464) 


(465) 


(466) 


(467) 


(468) 
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Then  substituting  equation  (468)  into  equation  (464), 


F„u  =  -(E|p'|2  +  com2)|p|"2x"'E/"2(|f|2  +  |pf)- 

•  [(m2-xE)2  +  E/2|p|2  sin2  0].  (469) 

Substituting  equations  (469),  (449),  (436),  and  (304)  into  equation  (462),  then 

F„,  =  -E/2(|f|2  +  |pf  )2(8|pf  |p|2  x2)-{4  -  12(E|p/|2  +  com2) 

•  [|p|2xE/2(|t|2  +  |^|2)]-‘[(m2-xE)2  +  E/2|p|2  sin2  0] 

-  4|p/|2(m2  -  xE)E''2(|t|2  +  |pf)-‘  +  2m2E/’2}.  (470) 

Simplifying  equation  (470),  then 

F„i  =  -(|t|2  +  rp'|2)(4rP'|2|p|V)-{|p|2x(2F2  +  m2) 

•  (|T|2  +  |V|2)  -  6(E|p'|2  +  com2)[(m2  -  xE)2 

+  E/2|p|2  sin2  6 ]  -  2|p/|2|p|2x(m2  -  xE)}.  (471) 

Using  equation  (43)  in  equation  (451), 

|t|2  +  |pf  =  |p|2  +  fpf  +  co2  -  2|p|co  cos  6 .  (472) 

Substituting  equation  (472)  into  equation  (471),  then 

F„i  =  -(|T|2  +  rP'|2)(4|J»'|2|p|4xT,{(2E'2  +  m2)x|p|2 

•  (|p|2  +  | P'12  +  co2  -  2|p|co  cos  6)  -  6(com2  +  EIp'I2) 

•  [(xE  -  m2)2  +  E/2|p|2  sin2  6}  +  2fpf|p|2x(xE  -  m2)} .  (473) 

Therefore,  substituting  equations  (448),  (449),  (460),  (473),  and  (436)  into  equation  (426),  then 
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Ia  =  E/L(2|p|  |p'|x)-,(m2E/"'x*‘  -  m2a>EE/",|p|-V  +  m2coE'''|p|-2 

+  x  'E2  (-2E'x  +  |p|2  +  |^|2  +  co2  -  2|p|  co  cos  0  +  2coE2  -  2|p|  E2  cos  0 
+  2E2  +  2E/2)  +  m2E2|p'|‘2E/~2|p|~4x‘3  { (2E/2  +  m2)x|p|2(|p|2  +  |^|2  +  a,2 
-  2|p|co  cos  9)  -  6(com2  +  E|pf  )[(Ex  -  m2)2  +  E'2|p|2  sin2  9] 

+  2Ip/|2|p|2x(xE  -  m2)}  -  (m2  -  xE)(2x3E2|p|2)-[8E2E2x  +  2co2|p|2  sin2  9 
+  2co2|p/|2  -  2co2E/2  -  2(E2  +  E22)(|p|2  +  |^|2  +  w2 . 2|p|a>  cos  9 
+  2coE2  -  2|p|E/  cos  0)]).  (474) 


Using  equation  (43), 


cos  0  =  (E  -  x)|p|-‘ 


(475) 


Substituting  equation  (475)  into  equation  (213),  then 


sin2  0  =  1  -  (E  -  x)2|p|-2. 


Then  substituting  equations  (475)  and  (476)  into  equation  (474), 


(476) 


la  = 


E/L(2|pJ  |p/|x)~1(m2E/“‘x~I  -  m2coEE/‘1|pr2x'1  +  m2a>E''‘|p|-2 
+  x  lE2  [-2E2x  +  |p|2  +  +  co2  -  2co(E  -  x)  +  2coE 2  -  2E'(E  -  x)  +  2E2  +  2E/2] 

+  m2E2 Ip2 |'2E/'2|p|'4x"3 { (2E/2  +  m2)x|p|2[|p|2  +  |p'|2  +  co2  -  2co(E  -  x)J 

-  6(com2  +  E|p/|2)[(Ex  -  m2)2  +  E'2|p|2  -  E22(E  -  x)2J  +  2|^|2|p|2x(xE  -  m2)} 

-  (m2  -  xE)(2x3E2|p|2)~‘{8E2E2x  +  2co2|p|2  -  2co2(E  -  x)2  +  2co2|p/|2  -  2co2E/2 

-  2(E2  +  E/2)[|p|2  +  |£f  +  a,2 . 2co(E  -  x)  +  2coE2  -  2E2(E  -  x)]}).  (477) 


Rewriting  equation  (477)  as  a  sum  of  powers  of  x~‘,  then 


la  =  E/L(2|p|  |p,|x)_,(K1  +  K2x-  +  K3x-2  +  K4X'3), 


(478) 
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where 


K,  =  m2aj(E/|p|2)-'  -  2  +  2a)E'-'  +  2  +  E(2E/|p|2)-(-2cu2),  (479) 

K2  =  m2E/'1  -  m2cuEE/"1|p|"2  +  E^^IpI2  +  |pf  +  co2  -  2coE  +  2coE'  -  2EE'  +  2E2  +  2E'2) 

+  2m2E2Q)|p,|'I|p|"4E/"2(2E/2  +  m2)|p|2  -  6m2E2fcr2|p|-4E'-2 

•  (cum2  +  E|p'|2)E2  +  6m2E2|p'|'2|p|'4E/_2(cum2  +  E|p'|2)E/2  +  2m2E3|p|‘2E/'2 

-  m2(2E/|p|2)'1(-2o)2)  +  E(2E/|p|2)*1[8E/E2  +  4co2E  -  2(E2  +  E/2)(2co  +  2E')],  (480) 

K3  -  m2E2|p/|‘2|p|'4E/'2(2E/2  +  m2)|p|2(|p|2  +  fpf  +  co2  -  2coE) 

+  12m2E2|pt2|p|-4E/-2(com2  +  E|pf  )m2E  +  m2E*fr|-*|p|-4E'_2 

•  [-6(com2  +  E|pf  )(2E/2E)  -  2|^|2|p|2m2]  -  m2(2E/|p|2)-[8E/E2  +  4co2E  -  4(E2  +  E/2)co 

-  4(E2  +  E/2)E2]  +  E(2E/|p|2)-,[2cu2|p|2  -  2cu2E2  +  2co2|p'|2  -  2co2E/2 

-  2(E2  +  E/2)(|p|2  +  |p'|2  +  co2  -  2coE  +  2o>E'  -  2EE/)],  (481) 

K4  =  -6m2E2|p,|'2|p|'4E/'2(com2  +  E|^|2)(m4  +  E'2fP|2  -  E'2E2) 

-  m2(2E/|p|2)-‘[2a)2|p|2  -  2cu2E2  +  2cu2|^|2  -  2co2E/2  -  2(E2  +  E'2) 

'  (|p|2  +  Ip'I2  +  co2  -  2coE  +  2cuE'  -  2EE2)]-  (482) 

Simplifying  equation  (479), 

Kj  =  a)E/"1|p|"2(m2  +  2|p|2  -  Ecu).  (483) 

Using  equations  (116)  and  (130)  in  equation  (483),  then 

Kt  =  a)E/-1|p|“2[m2  +  2E2  -  2m2  -  E(E  -  E')]  (484) 
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or 


K,  =  wE/_l|p|-2(E2  +  EE'-m2). 

Further,  simplifying  equation  (480), 

^2  =  K2l  +  K22  +  K23  +  K24, 

where 

K2,  =  2m2E2o>(2E/2  +  m2)|^|-2|p|-2E/-2  -  6m2E4(wm2  +  E|^|2)|^|-2|p|-"E/'2 
+  6m2E2|p/|~2|p|'4(a>m2  +  E|^|2)  +  2m2E3|p|-2E'~2, 

K22  =  |p|"2(4E3  -  2E3  -  2EE/a>  -  2EE/2), 

K2j  =  E/',|p|'2[m2|p|2  -  m2a>E  +  m2a>2  +  2a>2E2  -  2E3a>], 

K24  =  E'  ‘(IpI2  +  |^|2  +  w2  -  2coE  +  2a>E'  -  2EE'  +  2E2  +  2E/2). 

Substituting  equations  (116),  (209),  and  (130)  into  equation  (490),  then 

K24  =  E/_1[E2  -  m2  +  E/2  -  m2  +  (E  -  E')2  -  2(E  -  E')E 
+  2(E  -  E')E'  -  2EE'  +  2E2  +  2E/2] 


(485) 


(486) 


(487) 

(488) 

(489) 

(490) 


(491) 


K24  =  2E'-‘(E2  +  E/2  -  m2). 


(492) 
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Substituting  equation  (130)  into  equation  (489), 

K22  =  |p|-2[2E3  -  2E(E  -  co)co  -  2E(E  -  co)2]  (493) 

or 

K22  =  2coE2|p|“2.  (494) 

Using  equation  (130)  in  equation  (489),  then 

K23  =  E/_1|p|"2{m2|p|2  +  co[-m2E  +  m2(E  -  E')  +  2E2(E  -  E')  -  2E3]}.  (495) 

Simplifying  equation  (495),  then 

K23  =  E/_1|p|_2[m2|p|2  -  o^ir^E'  +  2E2E/)]  (496) 

or 


K23  -  m2E/_l|p|"2(|p|2  -  coE')  -  2coE2|p|~2. 


(497) 


Next,  simplifying  equation  (487), 


K21  =  2|p/|*2|p|"4E/"2[m2E2|p|2a>(2E/2  +  m2)  -  3m2E4(com2  +  E|J/|2) 

+  3m2E/2E2(com2  +  E|pf )  +  m2|^|2  |p|2E3  ],  (498) 


or  rearranging  terms,  then 


K21  =  2m2E2|p/|"2|p|“4E/-2[2cL)(|p|2  +  m2)E/2  +  m2coE/2  +  m2co|p|2  -  3m2coE2 

-3|^|2E3  +  3  |p/|2EE/2  +  |p/|2|p|2E].  (499) 

Substituting  equation  (116)  in  the  first  term,  then  equation  (499)  becomes 
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K21  =  2m2E2|p/|~2|p|~4E2~2[2o>E2E22  +  m2wE22  +  m2w|p|2  -  3m2wE2  -  3|^|2E3 

+  3|p/|2EE/2  +  |p/|2|p|2E],  (500) 

Substituting  equations  (130),  (116),  and  (209)  into  equation  (500),  then 

K21  =  2m2E2|^|-2|p|-"E/’2[2(E  -  E/)E2E/2  +  m2E'2(E  -  E2)  +  m2(E  -  E2)(E2  -  m2) 

-  3m2E2(E  -  E2)  -  3E3(E/2  -  m2)  +  3EE/2(E/2  -  m2)  +  (E/2  -  m2)(E2  -  m2)E],  (501) 

Simplifying  equation  (501), 

K21  =  2m2E2|p/|‘2|p|‘4E2'2(-2E2E/3  -  3m2EE/2  -  m2E/3  +  2m2E2E2  +  m4E2  +  3EE/4).  (502) 

Combining  the  first  and  fourth,  second  and  sixth,  and  third  and  fifth  terms,  then 

K21  =  2m2E2|p/|"2|p|'4E2'2[2E2E2(-E/2  +  m2)  +  3EE/2(-m2  +  E/2)  +  m2E2(m2  -  E/2)].  (503) 

Using  equation  (209)  in  equation  (503),  then 

K21  =  2m2E2|p/|'2|p|'4E2‘2(-2E2E2|p/|2  +  3EE22|pf  -  m2E2|^|2)  (504) 

or 

K2,  =  2m2E2|p|"4E2''(-2E2  +  3EE2  -  m2).  (505) 

Finally,  substituting  equations  (505),  (494),  (497),  and  (492)  into  equation  (486),  then 

K2  =  2m2E2|p|"4E2‘1(-2E2  +  3EE2  -  m2)  +  2wE2|p|-2  +  m2E2-'|p|-2(|p|2  -  wE2) 

-  2wE2fP|-2  +  2E2“'(E2  +  E22  -  m2)  (506) 
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or 


K2  =  m2E/_‘|p|‘2(|p|2  -  cuE')  +  2E/-'(E2  +  E'2  -  m2)  +  2m2E|p|-4E'-lF*2,  (507) 


where 


F*2  =  (-2E2  +  3EE2  -  m2)E. 


(508) 


Substituting  equation  (130)  into  equation  (508),  then 


F„  =  (E2  -  3 Ecu  -  m2)E. 


(509) 


Using  equation  (116),  equation  (509)  becomes 


Fjt2  =  (|p|2  -  3Ecu)E. 


(510) 


Using  equation  (116),  equation  (510)  becomes 


F*2  =  |p|2E  -  3(|p|2  +  m2)cu 


(511) 


or 


Fjn  =  -(3com2  -  |p|2E)  -  3|p|2cu. 


(512) 


Substituting  equation  (512)  into  equation  (507),  then 


K2  =  -2m2E|p|"4E/”,(3cum2  -  |p|2E)  +  f^E7-1  |p|-2. 


(513) 


where 


f*r2  =  2m2E|p|‘2(-3|p|2cu)  +  m2(|p|2  -  coE')  +  2(E2  +  E/2  -  m2)|p|2 


(514) 
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or 


f*2  =  -6m2Ea;  -  m2|p|2  -  ndcoE'  4-  2(E2  4-  E/2)|p|2.  (515) 

Using  equations  (130)  and  (116)  in  equation  (515),  then 

=  -6m2E(E  -  E')  -  m2(E2  -  m2)  -  m2(E  -  E')E'  +  2(E2  +  E/2)(E2  -  m2)  (516) 

or 

f„  =  -6m2E2  +  6m2EE/  -  m2E2  4-  m4  -  m2EE/  +  m2E/2  +  2E4  -  2m2E2  +  2E2E/2  -  2m2E/2  (517) 


or 


f„  =  -9m2E2  +  5m2EE/  +  m4  -  m2E/2  +  2E4  +  2E2E/2 

=  2E2(E2  +  E/2)  +  m4  -  m2(9E2  -  5EE2  +  E/2) 

=  2E2(E2  +  E/2)  +  m4  -  m2(7E2  +  2Eco  -  3EE2  +  E/2) 

=  2E2(E2  +  E/2)  -  m2(7E2  -  3EF  +  E/2)  +  m4  -  2m2Eo;.  (518) 


Substituting  equation  (518)  into  equation  (513),  then 

K2  =  -2m2Eo)E/"‘|p|"2  -  2m2E|p|“4E/“'(3com2  -  |p|2E) 

+  E/-,|p|-2[2E2(E2  +  E/2)  -  m2(7E2  -  3EE'  +  E/2)  +  m4].  (519) 

Combining  the  first  two  terms  of  equation  (519),  then 

K2  =  -2m2E|p|-4E'-'[3a)m2  -  |p|2(E  -  to)]  +  E/-*|p|-2[2E2(E2  +  E/2) 

-  m2(7E2  -  3EE2  +  E/2)  +  m4].  (520) 
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Using  equation  (130)  in  the  second  term,  then 


K2  =  -2m2E|p|~4E/~l(3cpm2  -  |p|2E')  +  E7"1  |p|-2[2E2(E2'  +  E/2) 

-  m2(7E2  -  3EE2  +  E/2)  +  m4].  (521) 

Equation  (481)  can  be  written  as 

K3  =  K3i  +  K32  +  K33  +  K34  +  K35,  (522) 

where 

K31  =  m2E2|p/|~2|p|'2E/'2(2E/2  +  m2)(|p|2  +  |^|2  +  gp2  -  2cpE),  (523) 

K32  =  12m4E3|p'|'2|p|'4E/'2(cpm2  +  Efpf),  (524) 

K33  =  m2E2|p/|"2  |p|'4E/‘2[-12(o)m2  +  E|^|2)E/2E  -  2|pf  IpP™2],  (525) 

K34  =  -m2(2E/|p|2)-’[8E/E2  +  4gp2E  -  4(E2  +  E/2)cp  -  4(E2  +  E/2)E2],  (526) 

K35  =  E(2E/|p|2)-,[2co2|p|2  -  2gp2E2  +  2gp2|^|2  -  2gp2E/2  -  2(E2  +  E/2) 

•  (IpI2  +  Ip'I2  +  cp2  -  2cpE  +  2gpE'  -  2EE')].  (527) 

Using  equation  (130)  in  equation  (523), 

K31  =  m2E2|p'|‘2|p|"2E/'2(2E/2  +  m2)[|p|2  +  +  (E  -  E2)2  -  2(E  -  E2)E]  (528) 

or 

K31  =  m2E2|p'|"2|p|"2E/‘2(2E/2  +  m2)(|p|2  +  |^|2  -  E2  +  E/2).  (529) 
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Substituting  equations  (116)  and  (209)  into  equation  (529),  then 


K31  =  2m2E2|p|~2E/~2(2E/2  +  m2). 


(530) 


Using  equation  (130)  in  equation  (524), 


K32  =  12m4E3|p'|'2|p|"4E/‘2[(E  -  E')m2  +  E|J/|2] 


(531) 


or 


K32  =  12m4E3|p/|'2|p|-4E/‘2[E(m2  +  |p'|2)  -  E'm2]. 


(532) 


Substituting  equation  (209)  into  equation  (532),  then 


K32  =  12m4E3|p'|"2|p|'4E/*2(EE/2  -  E'm2) 


(533) 


or 


K32  =  12m4E3|p/|"2|p|"4E/"'(EE/  -  m2). 


(534) 


Using  equation  (130)  in  equation  (525), 


K33  =  m2E2|p/|‘2|p|‘4E/'2{-6[(E  -  E')m2  +  E|pf]  (2E'2E)  -  2|^|2|p|2m2 } ,  (535) 


or 


K33  =  m2E2|p/|"2|p|"4E/*2[-12E/2E(Em2  -  E'm2  +  E|^|2)  -  2|^|2|p|2m2].  (536) 


Using  equation  (209),  equation  (536)  becomes 


K33  =  m2E2|p,|'2|p|‘4E/'2[-12E/3E(EE/  -  m2)  -  2m2|^|2|p|2]  (537) 
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or 


K33  =  12m2E3E/|p'|“2|p|"4(m2  -  EE7)  -  2m4E2|p|-2E7-2.  (538) 


Using  equation  (130)  in  equation  (526), 

K34  =  -m2(2E/|p|2)-‘[8E/E2  +  4(E  -  E')2E  -  4(E2  +  E/2)(E  -  E')  -  4(E2  +  E/2)E']  (539) 


or 

K34  =  -m2(2E/|p|2)-'(8E/E2  +  4E3  -  8E2E2  +  4E/2E  -  4E3  -  4E/2E),  (540) 

or  combining  terms. 


K34  —  0. 


(541) 


Further  using  equations  (116),  (209),  and  (130)  in  equation  (527), 

K35  =  E(2E/|p|2)-‘[-4cu2m2  -  2(E2  +  E'2)(|p|2  +  |^|2  -  co2  -  2EE2)].  (542) 


Again  using  equations  (130),  (209),  and  (116)  in  equation  (542), 

K35  =  E(2E7|p|2)_1[-4m2(E2  -  2EE7  +  E/2)  -  2(E2  +  E72)(-2m2)]  (543) 

or 


K35  =  4m2E2|p|*2. 


(544) 


Then  substituting  equations  (530),  (534),  (538),  (541),  and  (544)  into  equation  (522), 
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K3  =  2m2E2|p|-2E'-2(2E/2  +  m2)  +  12m4E3|^|-2|p|-',E/-(EE'  -  m2) 

+  12m2E3E/|p/|‘2|p|~4(m2  -  EE2)  -  2m4E2|p|-2E/-2  +  4m2E2|p|-2, 


(545) 


or  combining  the  second  and  third  terms  in  equation  (545)  and  using  equation  (209), 


K3  =  EE/-,|p|*2[2m2EE/'1(2E/2  +  m2)  -  12m2E2|p|-2(EE'  -  m2)  -  2m4EE/"  +  4m2EE/]  (546) 


or 


K3  =  2m2EE/'1  |p|'2[ EE/'I(2E/2  +  m2)  -  6E2|p|-2(EE'  -  m2)  -  m2EE/'1  +  2EE2]  (547) 


or 


K3  =  2m2EE/-‘|p|-2[4EE/  -  6E2(EE2  -  m2)|p|-2]. 


(548) 


Equivalently, 


K3  =  4E2m2E,-‘|p|-4[2E/|p|2  -  3E(EE2  -  m2)]. 


(549) 


Using  equations  (130)  and  (116)  in  equation  (549),  then 


K3  =  4E2m2E/-‘|p|-4[2E/|p|2  +  3m2(E2  +  <o)  -  3(|p|2  +  m2)E2]  (550) 


or 


K3  =  4E2m2E/-1|p|-4(3cum2  -  E2|p|2). 


(551) 


Multiplying  terms  in  equation  (482),  then 


K,  =  K*.  +  K,2, 


(552) 
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where 


K,i  =  -6m2E2|p'|'2|p|~4E2'2(a>m6  +  com2E/2|p|2  -  com2E/2E2  +  m4E|J/|2  - 

+  EE/2^|2|p|2  -  E3E/2|p/|2,  (553) 

K42  =  -mVE'“  +  m2co2E2E/“|p|-2-m2co2|^|2E/"|p|-2  +  m2co2E2fP|-2  +  m’E’E2'1 
+  m2E2|p'|2E/'1|p|~2  +  E2m2a>2E2"‘|p|"2  -  2m2a>E3E2~'|p|"2  +  2m2coE2|p|-J 
-  2m2E3|p|‘2  +  m2E/  +  m2E2|p/|2|p|"2  +  E/m2co2|p|"2  -  2m2coEE/|p|’2 

+  2m2cpE/2|p|"2  -  2m2E/2E|p|"2.  (554) 

Further,  by  multiplying  in  equation  (553),  then 

K41  =  -6com8E2|p,|‘2|p|"4E/'2  -  6m4E2co|p|2|^r2|p|‘4 

+  6m4E4co|p/|"2  |p|'4  -  6m6E3|p|'4E2~2  -  6m2E3|p|~2  +  6m2E8|p|"4.  (555) 

Using  equation  (116)  to  reexpress  factors  of  E2  in  the  third  and  sixth  terms  of  equation  (555),  then 

K41  =  -6com8E2|p'|~2|p|~4E/"2  -  bm4E2co|p|2|^|-2|p|-4  +  6m4E2 
•  (|p|2  +  m2)co|p'|‘2|p|"4  -  6m6E3|p|"4E/"2  -  6m2E3|p|-2 
+  6m2E3(|p|2  +  m2)|p|'4.  (556) 

Simplifying  and  using  equation  (209),  equation  (556)  becomes 

K,,  =  -6com8E2|p/ 1-2 |p|_4E/_2  +  6m6E2a>|p'|~2|p|~4  -  6m4 

•  (E22  -  |p,|2)E3|p|'4E2'2  +  6m4E3|p|-4  (557) 

or 

K41  =  -6m4E2|p'|"2|p|"4E2"2(com4  -  m2coE22  -  E|J/|4).  (558) 


89 


Using  equation  (209)  in  the  second  term  of  equation  (558),  then 


K4,  —  6m4E2|p|~‘,E2~2(a;m2  +  E|p'|2). 


(559) 


Next,  combining  terms,  equation  (554)  becomes 


K« 


-m2(E  -  E2)2E2"‘  +  2m2co2E2E/'1|p|"2  -  m2co2|p'|2E2'1|p|-2  +  2m2w2E/|p|'2  +  m2E2E/" 


+  m2(E2  +  oj)2|p/|2E/'l|p|"2  -  2m2coE2(E/  +  co)E/_1|p|-2  +  2m2coE2|p|-2  -  2m2E3|p|'2 
+  m2E2  +  m2E/|p/|2|p|~2  -  2m2co(E2  +  co)E2|p|-2  +  2m2wE/2|p|-2  -  2m2E/2E|p|-2,  (560) 


where  equation  (130)  has  been  used  in  the  1st,  6th,  7th,  and  12th  terms.  Further,  combining  terms  in 
equation  (560),  then 


K42  =  2m2E  -  m2co2|p,|2E/~1|p|~2  +  m2E2|p/|2E/'1|p|-2  -  2m2E3|p|* 


+  m2E/|p/|2|p|"2  -  2m2E/2E|p|"2 


(561) 


or 


K«  =  m2(2E  -  co2|p'|2E/‘1|p|~2  +  E2|p'|2E/- |p|"2  -  2E3|p|"2 

+  EV|2|p|-2  -  2E/2E|p|"2).  (562) 

Using  equation  (130)  and  substituting  for  co  in  the  second  term  of  equation  (562),  then 

K42  =  2m2E|p|-2(|p|2  +  fcf  -  E2  -  E22).  (563) 

Substituting  equations  (116)  and  (209)  into  equation  (563),  then 

K42  =  -4nVE|p|~2.  (564) 
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Finally,  substituting  equations  (559)  and  (564)  into  equation  (552),  then 

K4  =  6m4E2|p|~4E2~2(com2  +  E|p'|2)  -  4m4E|p|"2  (565) 

or 

K4  =  2m4E|p|‘4E2"2(3Ecom2  +  3E2|p'|2  -  2|p|2E'2).  (566) 

Substituting  equation  (209)  into  the  second  term  of  equation  (566)  and  using  equation  (130),  then 
K4  =  2m4E|p|~4E2"2[3Ea>m2  +  3E2E/2  -  3E(E2  +  co)m2  -  2|p|2E/2]  (567) 

or 

K4  =  2m4E|p|-4E/-'(-3Em2  +  3E2E2  -  2|p|2E2).  (568) 

Using  equation  (130)  in  the  first  term  of  equation  (568),  then 

K4  =  2m4E|p|-4E/-,[-3com2  +  3E'(E2  -  m2)  -  2|p|2E2].  (569) 

Finally,  substituting  equation  (116)  into  equation  (569),  then 

K4  =  -2m4E|p|'4E/"1(3ojm2  -  |p|2E/).  (570) 

Then  by  substituting  equations  (485),  (521),  (551),  and  (570)  into  equation  (478), 

I.  =  E/L(2|p|  |p'|x)',(ajE/‘1|p|'2(E2  +  EE2  -  m2)  +  x"  {-2m2E|p|-4E/-' 

•  (3com2  -  |p|2E0  +  E/_* |p|_2[ 2E2(E2  +  E22)  -  m2(7E2  -  3EE2  +  E22)  +  m4]} 

+  x-2[4E2m2E2-'|p|-4(3com2  -  E2|p|2)]  +  x-2[-2m4E|p|-4E2- 

•  (3com2  -  |p|2E2)]).  (571) 
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Collecting  terms  with  the  common  factor  (3com2  -  |p|2E7),  equation  (571)  can  be  rewritten  as 


I.  =  L(2|f»|  |p'|)"{(3wm2  -  |p|2E2)|p|-4x-4Fo1  +  |p|-Vw(E2  +  EE2  -  m2) 

+  |p|-V2[2E2(E2  +  E/2)  -  m2(7E2  -  3EE2  +  E/2)  +  m4]},  (572) 

where 

Fal  =  -2m2Ex2  +  4E2m2x  -  2m4E.  (573) 

Equation  (573)  may  be  rewritten  as 

Fal  =  -2m2E(>c2  -  2Ex  +  m2) 

or 

F„,  =  -2m2E[(E  -  k)2  -  E2  +  m2].  (574) 

Using  equations  (43)  and  (116)  in  equation  (574),  then 

Fal  =  -2m2E(|p|2  cos2  6  -  |p|2).  (575) 

Using  equation  (213),  equation  (575)  becomes 

Fal  =  2m2E|p|2  sin2  0.  (576) 

Substituting  equation  (576)  into  equation  (572),  then 

I„  =  L(4|p|  |p'|)'l{4Em2  sin2  0(3com2  -  |p|2E')|p|-2>r4  +  [4E2(E2  +  E/2)  -  2m2 

•  (7E2  -  3EE'  +  E/2)  +  2m4]  rP|-2>c-2  +  2co(E2  +  EE2  -  m2)|p|-V‘}.  (577) 
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Proceeding  also  to  reduce  equation  (425),  one  first  rewrites  equation  (425)  in  the  form 

I*  =  Im  +  1*2  +  1*3  +  1*4  +  1*5  +  Its,  (578) 

where 

I*i  =  FtiU,,  (579) 

F»i  =  -[2 x~'(cu  -  |p|  cos  6)  4-  |p|2k"2  sin2  8  -  l]|p/|'1|a|'1 

•  (|T|2  +  Ip'I2)'1  -  2cu2E/X-irpt,|a|-' 

•  (|T|2  +  |p'|2)-2a-6|a|-2,  (580) 

Lt,  =  In  [(1  4-  Ip'I  |a|)(l  -  |p'|  |a|)->],  (581) 

It2  =  -*‘,E/'l|p'|‘,|6|"1  in  1(1  4-  Ip'I  |6|)(1  -  Ip'I  |6|)_1],  (582) 

Ita  =  m2(|T|2  +  |p'|2)-'E'-2(4BA-‘D-2),  (583) 

It4  =  Ib41  +  It42  4"  It43  4"  It44,  (584) 

Im,  =  -4m2E2E/"2(|f|2  4-  |p/|2)“2(A  4-  B  4-  C)-2[8D-2(1  4-  B/3C)(2A  4-  B) 

+  16/3C  4-  (8A/3C)(2C  4-  B)D"2],  (585) 

I»42  =  -4m2E2E'“2(|f|2  4-  |p/|2)“2(A  4-  B  4-  C)'1  [24D-4A-'(1  4-  B/3C)B 

•  (D2  4-  2AC  4-  BC)  4-’l6D-4(D2  4-  2AC  +  BC)],  (586) 

1*43  =  -4m2E2E,_2(|f|2  4-  |p'|2)-2[-16D-2A-'(l  +  B/3C)  -  8B(3CD2A)-  +  8D'2],  (587) 
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U.  =  -4(B  +  2C)x',E/'1(|T|2  +  |pf  )-2D-2(A  +  B  +  C)-'[8E'E2x  +  2co2|p|2  sin2  0 
+  2a>2|p,|2  -  2a>2E/2  -  2(E2  +  E/2)(|p|2  +  |pf  +  co2  -  2|*p|a;  cos  0 

+  2cjE2  -  2|p|E2  cos  0)],  (588) 

lbs  =  2(|T|2  +  |pf  )-2(l  -  |^|2  |a|2)-[4x-(E2  +  E/2)(co  -  |p|  cos  0)  -  4E2 

+  4|p|2E/2x'2  sin2  0  +  4E/co2x-‘],  (589) 

U  =  -4a>2E/x‘1(|f|2  +  |^|2)-2(1  -  a  ■  fe|a|-2)(l  -  ff/|2  |a|2)-\  (590) 

Substituting  equation  (244)  into  equation  (581), 

Lm  =  in{[l  +  2|p'|  |T|(|T|2  +  |^|2)-‘][1  -  2|p'|  |T| 

•  (|T|2  +  Ip'I2)-]'’ }  (591) 

or 

Lm  =  2  in  (<|T|  +  |p'|)(|T|  -  |p'|)-'].  (592) 

To  reduce  equation  (580),  one  first  notes  that,  using  equations  (244),  (245),  and  (41), 

a •  S  | a |-2  =  T-  C(|T|2  +  |pf  )(2coE/|t|2)-*.  (593) 

Substituting  equation  (300)  into  equation  (593),  then 

a  •  6  | a | ~2  =  (4coE/|f|2)-' 

•  (|pf  +  2FCO-  |T|2)(|T|2  +  Ip'I2).  (594) 

Substituting  equations  (244)  and  (594)  into  equation  (580),  then 

Fm  =  F.„  +  Ffi  12/  (595) 
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where 


Fmi  =  -[2x~'(w  -  |p|  cos  6)  +  |p|2x  2  sin2  9  -  l]|p'|'1|a|"1 

•  (|T|2  +  l^l2)-1  +  co(4x|H  |t|)-,  (596) 

Fmz  =  -co(|p'|2  +  2E/co)(4x|p/|  |f|3)-'.  (59 7) 

Substituting  equations  (209)  and  (130)  into  equation  (597), 

Fm  =  -co(E/2  -  m2  +  2Eoj  -  2co2)(4x|p/|  |T|3)-1.  (598) 

Substituting  equation  (130)  into  equation  (598),  then 

Fm2  =  -co(E2  -  m2  -  oj2)(4x|p/|  |T|3)'\  (599) 

or  again  using  equation  (116), 

FM2  ”  -co(|p|2  -  co2)(4x|p/|  |T|3)'‘.  (600) 

Further,  using  equations  (43)  and  (244)  in  equation  (596),  then 

Fm  =  -{2x-‘(co  -  E  +  x)  +  |p|V2[l  -  (E  -  x)2|p|-2]  - 1}(2|^|  |t|)- 

+  o)(4x|p/|  |T|)-'  (601) 

or 

FM1  =  -(2|^|  |t|)'‘[|p|V2  -  E2x"2  +  3oj(2x)“'],  (602) 

or  using  equation  (116), 

Fm  =  -(2|p/|  |t|)‘1[3oj(2x)‘1  -  m2x‘2].  (603) 
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Substituting  equations  (603)  and  (600)  into  equation  (595),  then 


Fm  =  -(2|p'|  |t|)_1[3o)(2k)"1  -  mV2]  -  co(|p|2  -  co2)(4x|J>/|  |f|3)_1 

or 

Fm  =  (8|p/|  |f|)_1[4m2x“2  -  6o)k~1  -  2o>(|p|2  -  co2)x“1|t|’2]. 
Then  substituting  equations  (605)  and  (592)  into  equation  (579), 

hi  =  (4|p/|  |T|)_1[4mV2  -  6cl)x_1  -  2o)x_1|f |“2(|p|2  -  co2)] 

•  M(|f|  +  Ip/I)(|t|  -  Ip'D"]. 

Next  substituting  equation  (245)  into  equation  (582),  then 

1.2  =  -x-‘E/-*E/|p'|"  in  [(1  +  |HE'-‘)(1  -  |p/|E/',)‘1] 


I42  =  -K"|P'I“  in  [(E'  +  |p'|)(E/  -  Ip'D-]. 

Next,  substituting  equations  (304),  (436),  and  (449)  into  equation  (583), 

Ii3  =  -m2|p|~2xf2  +  E|p|'2x—. 

Equation  (585)  may  be  rewritten  as 

Imi  =  -4m2E2E/-(|T|2  +  |p,|2)"2[8(B  +  3C)(2A  +  B)  +  16D2 
+  8A(2C  +  B)][3CD2(A  +  B  +  C)2]- 


or 


(604) 


(605) 


(606) 


(607) 


(608) 


(609) 


(610) 
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IM1  =  -4m2E2E/_2(|T|2  +  fP'|2)-'(24AB  +  8B2  +  64AC 
+  24BC  +  16D2)[3CD2(A  +  B  +  C)2]-. 


(611) 


Substituting  equation  (283)  in  the  numerator  of  equation  (611),  then 

IM.  =  -32m2E2E'-2(|T|2  +  |pf  )-2BC-D-2(A  +  B  +  C)-.  (612) 

Next,  simplifying  equation  (587), 

U  =  -4m2E2E/"2(|f  |2  +  |pf  )-2Ffc43,  (613) 

where 

F*43  =  -16(1  +  B/3C)D_2A_1  -  8B(3CD2A)-  +  8D~2  (614) 

or 

Fm3  =  8D~2A_1(-2  -  BC"1  +  A).  (615) 

Substituting  equation  (615)  into  equation  (613),  then 

1*43  =  -32m2E2E/-2(|T|2  +  |p'|2)~2D~2A",(-2  -  BC^1  +  A).  (616) 

Equation  (588)  may  be  rewritten  as 

lb44  =  -4(B  +  2C)x_1E/'1(|f|2  +  |pf)-2D-2(A  +  B  +  C)-’FM4,  (617) 

where 

F644  =  8E7E2x  +  2o>2|p|2  sin2  0  +  2co2|pf  -  2w2E'2  -  2(E2  +  E/2)(|p|2  +  |pf  +  o>2 

-  2|p|ca  cos  0  +  2coE/  -  2|p|E7  cos  0).  (618) 
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Using  equations  (130),  (43),  (116),  (213),  and  (209)  in  equation  (618),  then 

FM4  =  SE'E2*  +  2(E  -  E')2(E2  -  m2)[l  -  (E  -  x)2(E2  -  m2)'1]  -  2(E  -  E2)2m2  -  2(E2  +  E/2)[E2  +  E/2 
-  2m2  +  E2  -  2EE2  +  E/2  -  2E(E  -  x)  +  2(E  -  EOE'],  (6: 

or  simplifying  equation  (619),  then 


Fk44  =  8m2EE2  -  2(E  -  E2)2*2.  (620) 

Substituting  equation  (620)  into  equation  (617),  then 


Im4  =  -4(B  +  2C)[8m2EE/  -  2(E  -  EOVHxE'dTl'2  +  |^|2)D2(A  +  B  +  C)]-‘. 


(621) 


Using  equations  (612)  and  (616),  then 


U  +  IM3  =  -32m2E2B[E/2(|T|2  +  |p'|2)2CD2(A  +  B  +  C)]"* 

-  32m2E2(-2  -  BC-  +  A)[E'2(|t|2  +  |^|2)2D2A]-  (622) 


or 


Im.  +  I„43  =  Im.3  -  32m2E2(-2  +  A)[E/2(|f|2  +  |pf  )2D2A]-‘, 

where 


Im,3  =  -32m2E2B[E/2(|T|2  +  |p'|2)2CD2]-'[(A  +  B  +  C)'1  -  A-'], 


Equation  (624)  may  be  rewritten  as 


Imu  =  32m2E2B(B  +  C)E/_2(|f|2  +  |pf  )'2 


•  [AD2C(A  +  B  +  C)]". 


(623) 


(624) 


(625) 


98 


Substituting  equation  (625)  into  equation  (623),  then 


U  +  Im3  =  32m2E2B(B  +  C)[E'2(|t|2  +  frf  )2AD2C(A  +  B  +  C)]" 

-  32m2E2(A  -  2)[E/2(|T|2  +  |^|2)2D2A]-'.  (626) 

Substituting  equations  (626)  and  (586)  into  equation  (584),  then 

IM  =  IMo  -  32m2E2(A  -  2)[E/2(|t|2  +  |pf  )*D*A]-,  (627) 

where 

Ifr4a  =  Ib4  b  +  h>44,  (628) 

IMfc  =  32m2E2B(B  +  C)[E/2(|t|2  +  |^|2)2AD2C(A  +  B  +  C)]“ 

-  32m2E2[E/2(|T|2  +  |^|2)2]-'{(1  +  B/3C)3B(D2  +  2AC  +  BC)[D4A(A  +  B  +  C)]"1 
+  2(D2  +  2AC  +  BC)[D4(A  +  B  +  C)]1}.  (629) 

Combining  terms  in  equation  (629),  then 

Ifc4b  =  32m2E2[E/2(|t|2  +  |p'|2)2ACD4(A  +  B  +  C)]“FMfc,  (630) 

where 

FMt  =  BD2(B  +  C)  -  3BC(1  +  B/3C)(D2  +  2AC  +  BC)  -  2AC(D2  +  2AC  +  BC).  (631) 
Equivalently,  equation  (631)  is 

FMt  =  B2D2  +  BCD2  -  (3BC  +  B2)(D2  +  2AC  +  BC)  -  2AC(D2  +  2AC  +  BC)  (632) 

or 
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Fm„  =  B2D2  +  BCD2  -  (D2  +  2AC  +  BC)(3BC  +  B2  +  2AC) 


(633) 


or 


or 


Fm*  =  D2(-2BC  -  2AC)  -  (2AC  +  BC)(3BC  +  B2  +  2AC) 


(634) 


Fm‘  “  -2CD2(A  +  B  +  C)  +  2CD2  -  (2AC  +  BC)(3BC  +  B2  +  2AC). 
Substituting  equation  (283)  into  the  second  term  of  equation  (635),  then 

Fm‘  =  -2CD2(A  +  B  +  C)  +  C(-B2C  -  8AC2  -  8ABC  -  2AB2  -  4A2C  -  B3), 


(635) 


(636) 


Fm6  -  -2CD2(A  +  B  +  C)  -  B2C(A  +  B  +  C)  -  ACD2  -  8AC2(A  +  B  +  C).  (637) 

Substituting  equation  (637)  into  equation  (630),  then 

U  =  32m2E2[E/2(|f|2  +  |p'|2)2ACD4(A  +  B  +  C)]-[-2CD2(A  +  B  +  C) 

-  B2C(A  +  B  +  C)  -  ACD2  -  8AC2(A  +  B  +  C)].  (638) 

Substituting  equations  (638)  and  (621)  into  equation  (628),  then 

L.,.  =  -64m2E2E/'2(|f|2  +  |f/|2)‘2A'*D'2  -  32m2E2B2E/'2 

'  (|T|2  +  |H2)-2A-D-4  -  256m2E2CE/'2(|f|2  +  |^|2)-2D-  -  Ib4c,  (639) 

where 

Imc  =  32m2E2E/  2(|f|2  +  |p/|2)-2D'2(A  +  B  +  C)‘‘  +  4(B  +  2C)[8m2EE/  -  2(E  -  E7)2x2] 

•  x-E'-‘(|T|2  +  |pf)‘2D-2(A  +  B  +  Q". 
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Rewriting  equation  (640),  then 


IMc  =  4{8xm2E2  +  E'(B  +  2C)[8m2EE/  -  2(E  -  E')2x2]} 

•  x'1E/~2[(|t|2  +  |f/|2)2D2(A  +  B  +  C)]-‘  (641) 

or 


IMC  =  4{8xm2E2  -  E'[(2A  +  B)  -  2(A  +  B  +  C)][8m2EE'  -  2(E  -  E')2x2]} 

•  x-E'-2(|T|2  +  |f/|2)'2D'2(A  +  B  +  C )-.  (642) 

Then  equation  (642)  becomes 

IMc  =  IM,  +  8E'[8m2EE'  -  2(E  -  E')2x2]x"E''2(|T|2  +  |J/|2)-2D-2,  (643) 

where 

\iAd  =  4{8xm2E2  -  E'(2A  +  B)[8m2EE'  -  2(E  -  E')2x2]} 

•  [xE/2(|t|2  +  |f/|2)2D2(A  +  B  +  C)]M.  (644) 


Using  equations  (276)  to  (278), 


A  +  B  +  C  =  1  -  fpf|6|2  -  2|^|2(a  •  £-|£|2) 
-  |pf  |a  -  6|2 


(645) 


or 


A  +  B  +  C  =  1  -  |p'|2|a|2. 

Substituting  equation  (244)  into  equation  (646),  then 

A  +  B  +  C  =  1  -  4fP'|2|T|2(|T|2  +  |H2)-2. 


(646) 


(647) 
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Simplifying  equation  (647),  then 


A  +  B  +  C  =  (|TM;>f)2(|T|2  +  |^|2)-2. 

Substituting  equation  (448)  into  equation  (648),  then 

A  +  B  +  C  =  (2cox)2(|t|2  +  |^|2)-2. 

Next,  substituting  equations  (304),  (436),  (449),  and  (649)  into  equation  (644),  then 

Im*  =  4{8xm2E2  -  E/[8m2EE/  -  2(E  -  E/)2x2][2m2E/~2  -  4|Jy|2(m2  -  xE)E/_2 

•  (|T|2  +  li^l2)'1]}  {xE/2(16|p'|2|p|2k2)[E/2(|t|2  +■  |p'|2)J]-‘4a)2x2}-1. 

Using  equation  (130)  and  simplifying  the  result,  then  equation  (650)  becomes 


(648) 


(649) 


(650) 


Imi  -  (|T|2  +  |p,|2)(4|p'|2|p|VE/ai2)'1{2xm2E/E2 

•  (|T|2  +  l^l2)  -  (4m2EE7  -  a;2x2)[m2(|f|2  +  |pf )  -  2|pf(m2  -  xE)]}.  (651) 


Using  equation  (448)  in  equation  (651),  then 


Imj  —  (cox  +  |p'|2)(2|p/|2|p|2x5E/co2)'IFMd, 


(652) 


where 


F64d  -  (2cox  +  2|p'|2)(2xm2E2E/  -  4m4EE/  +  m2co2x2)  +  2|pf  (m2  -  xE)(4m2EE/  -  co2x2). 


(653) 


Multiplying  in  equation  (653),  then 


4cox2m2E2E/  -  8com4EE'x  +  2m2co3x3  +  4|p'|2m2E2E/x  -  8m4|^|2EE/  +  2m2|p'|2co2x2 


+  8m4EE/|p/|2  -  2|pf  m2co2x2  -  8m2|J/|2E2E/x  +  2|pf  co2Ex3 


(654) 
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or 


F„d  =  4a>m2E2E/Jc2  +  2m2w¥  +  2fp'|2a>2E>c3  -  8wm4EE'>c  -  4|p'|2m2E2E'>c. 
Substituting  equation  (655)  into  equation  (652),  then 

I64d  =  (cox  +  |p'|2)(2|p'|2|p|2JC5E/a)2)‘,(4a)m2E2E/JC2  +  2m2a)V  +  2fP'|2«2E>c3 

-  8wm4EE/jc  -  4|p,|2m2E2E'jc). 

Substituting  equations  (652)  and  (655)  into  equation  (643),  then 

IMc  =  (cox  +  |pf)(2|pf  |p|VE/w2)-(4wm2E2E')c2  +  2m2wV  +  2|pf  co2Ex3 

-  8wm4EE/)c  -  4|p'|2m2E2E/)c)  +  8E'[8m2EE'  -  2(E  -  E')2*2] 

•  [jcE/2(  |T|2  +  |p'|2)2D2]'1. 

Substituting  equation  (657)  into  equation  (639),  then 

u  =  -64m2E2[E/2(|T|2  +  |p/|2)2AD2]'1  -  32m2E2B2 

•  [E/2(|t|2  +  |p,|2)2AD4]"1  -  256m2E2C[E/2(|T|2+  |^|2)2D4]-‘ - (cojc  +  fc|2) 

•  (2|p,|2|p|2K5E/w2)'1(4wm2E2E/jc2  +  2m2wV  +  2|pf  «2Ex3 

-  8com4EE/)c  -  4fP'|2m2E2E'>c)  -  8E/[8m2EE/  -  2(E  -  E')V] 

•  [>cE/2(|t|2  +  |p'|2)2D2]'1. 

Substituting  equation  (658)  into  equation  (627),  then 


(655) 


(656) 


(657) 


(658) 
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Im  =  -64m2E2[E/2(|f|2  +  |p'|2)2AD2]-‘  -  32m2E2B2 

•  [E/2(|T|2  +  |pf  )2AD4]-  -  256m2E2C[E/2(|T|2  +  )2D4]- 

-  8E/[8m2EE/  -  2(E  -  E')2x2][xE'2(|t|2  +  |pf)2D2]-'  -  (cox  +  |pf ) 

'  (2|p'|2|p|2x5E/co2)'1(4com2E2E/x2  +  2m2co3x3  +  2|p,|2co2Ex3 

-  8com4EE/x  -  4|^|2m2E2E'x)  -  32m2E2(A  -  2) 

•  [E'2(|T|2  +  |pf)2D2A]-‘.  (659) 

To  reduce  equation  (589),  one  first  notes  using  equation  (244)  that 

1  -  |p'|2|a|2  =  1  -  (2|p'|  |t|)2(|t|2  +  Ip'I2)-2,  (660) 

or  substituting  equation  (448),  then 

1  -  |p,|2|a|2  =  (2cox)2(|T|2  +  |pf  )-2.  (661) 

Substituting  equations  (661)  and  (43)  into  equation  (589),  then 

Im  =  2(2cox)-2{4EW  +  4E/2cox-  -  4(E2  +  E/2)x-(E  -  x)  -  4E2  +  4|p|2E/2x-2 

•  [l  -  (E  -  x)2|pp2]  +  4E/co2x'1} .  (662) 

Simplifying  equation  (662),  then 

Im  =  2co-2x-2(E2cox"  +  E/2cox-‘  -  E3x“  -  E/2Ex-  +  |p|2E/2x“2  -  E/2E2x'2 

+  2E/2Ex"‘  +  E/co2x"1).  (663) 

Substituting  equation  (116)  into  equation  (663),  then 

I.s  =  2co"2x"3(E2co  +  E/2co  +  E'co2  -  E3  +  E/2E  -  nTE'V).  (664) 
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Using  equation  (130)  in  equation  (664),  then 

I„  =  -2m2E/2cu-2x-4.  (665) 

Next  substituting  equations  (594)  and  (660)  into  equation  (590),  then 

I66  =  -4cu2E/x-‘(|T|2  +  |pfHl  -  (4o)E/|t|2)"i(|p/|2  +  2E/cu 
-  |T|2)(|T|2  +  |pf)][l  -  4|p,|2|T|2 

’  (|T|2  +  Ip'I2)'2]'1  (666) 

or 

U  =  -cu[x(|T|2  -  Ppf  )2|T|2]-‘(-|pf  +  2E'cu|T|2 

-  2E/copPf  +  |T|4).  (667) 


Substituting  equation  (448)  into  equation  (667),  then 

IM  =  [|pf  (2x2|T|2)-‘(ppf  +  2E/co)  -  EW2  -  |T|2(2x2)-1  ](2cox)  2 .  (668) 

Using  equation  (451)  in  equation  (668),  then 

I4t  =  (2cox)"l[|p'|2(2x2|t|2)",(|f/|2  +  2E/co)  -  EW2  -  [P^2  -  2co(E  -  x)  +  co2](2x2)-].  (669) 

Substituting  equation  (130)  into  the  second  term  of  equation  (669),  then 

Ii6  =  (2cox)'I[|p'|2(2x2|T|2)'l(|p/|2  +  2E/co)  -  cox'1  -  (|p|2  -  co2)(2x2)'‘]  (670) 


or 


I66  =  I„61  -  (2x2)-‘  -  (|p|2  -  co2)(4cox3)-1, 


(671) 
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where 


1*6,  =  |^|2(|^|2  +  2E/o))(4a)K3|f|2)"1.  (672) 

Using  equations  (448),  (209),  (130),  and  (116)  in  equation  (672),  then 

I*6i  =  (|T|2  -  2o)x)(|p|2  -  o)2)(4o)x3|f I2)'1  (673) 

or 

I*6i  =  -(|p|2  -  o)2)(2jc2|t|2)"‘  +  (|p|2  -  a)2)(4a)x3)''.  (674) 

Substituting  equation  (674)  into  equation  (671),  then 

1*6  =  -(|p|2  -  o)2)(2x2|f|2)-'  -  (2k2)"'.  (675) 

Then  using  equations  (609),  (659),  (665),  and  (675),  one  obtains 

1*3  +  1*4  +  1*5  +  1*6  =  1*36,  (676) 

where 

1*36  —  _(|p|2  -  co2)(2jc2|T|2)  1  -  (2jc2)"'  -  m2|p|~2x:~2  +  EIpI'2*:-1  -  2m2E/2ar2x:'4 

+  1*361  +  1*362  +  1*363  +  1*364  +  1*365  +  1*366,  (677) 

1*36.  =  -32m2E2(A  -  2)[E/2(|T|2  +  |^|2)2D2A]-,  (678) 

1*362  =  -64m2E2[E/2(|T|2  +  |p,|2)2AD2]'1,  (679) 

1*363  =  -32m2E2B2[E/2(|f|2  +  |^|2)2AD4]-,  (680) 

1*364  =  -256m2E2C[E/2(|T|2  +  |^|2)2D4]-,  (681) 
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16365  =  -8E/(8m2EE/  -  2a,V)[xE'2(|T|2  +  |p'|2 )2D2]-‘,  (682) 

where  equation  (130)  has  been  used,  and 

16366  =  -(2cjx  +  2|p'|2)(4|p/|2|p|2x5E/co2)'1(4com2E2E/x2  +  2m2co3x3b  2|pfco2Ex3 

-  8wm4EE/x  -  4|p'|2m2E2E/x).  (683) 

Substituting  equations  (304)  and  (449)  into  equation  (678),  then 

I6361  =  -32m2E2(m2E/"2  -  2)E2-2[l6|p'|2E/-2|p|2x2m2E/-2]-'  (684) 

or 

I6361  =  -2E2m2[|pf  |p|2x2]-‘  +  4E2E/2(|pf  |p|2x2)-.  (685) 

Substituting  equations  (304)  and  (449)  into  equation  (679), 

It362  =  -4E2E/2(|p/|2|p|2x2)''.  (686) 

Substituting  equations  (304),  (449),  and  (436)  into  equation  (680), 

Im«3  =  -2m4E2(|p|4x4)'1  +  4m2E3(|p|4x3)-’  -  2E4(|p|4x2)-.  (687) 

Substituting  equation  (446)  and  (449)  into  equation  (681)  and  using  equations  (43)  and  (213),  then 

It364  =  4m6E2(|p/|2|p|4x4)'1  -  8m4E3(|pf  |p|4x3)-  +  4m2E4 

•  (|pf|p|V)"  +  4m2E2E/2(|p'|2|p|2x4)'1  -  4m2E/2E4 

•  (|p,|2|p|4x4)"*  +  8m2E3E/2(|p/|2|p|4x3)'1  -  4m2E2E/2 

•  (IpflpIV)-.  (688) 
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Next,  substituting  equation  (449)  into  equation  (682), 


1*365  =  -4m2EE/2(|p'|2|p|V)-1  +  E/o)2(|pf  |p|2,c)- 


(689) 


Equation  (683)  is 


1*366  =  -2m2E2(|pf  |p|V)-  -  m2w2(|p,|2|p|2E/x)'1  -  Ea>(|ppE'ic)-‘ 

+  4m4E(|p'|2|p|V)-.  +  2m2E2(|p| W)-  -  2m2E2(|p|  W)- 
-  m2o)(|p|2E2x2)-‘  -  |pf  E(|p|2EV)"  +  4m4E(|p|W)- 
+  2|p,|2m2E2(|p|2o)2x4)"1. 


Substituting  equations  (685)  to  (690)  in  equation  (677)  and  collecting 
powers  of  x,  then  6 


(690) 

terms  containing  various 


1*36  -  Ifc.iJC'1  +  lbx2x~2  +  \bK3x~3  +  lbK4x~\ 

where 

I*x.  =  E|p|-2  -  m2o)2(|p'|2|p|2E/)'1  -  EwdppE2)-* 

+  E/o)2(|p/|2|p|2)'1, 

I*,2  =  -|p|2(2|f|2)-‘  +  o)2(2|T|2)-‘  -  V2  -  m2|p|~2  -  4E2m2 
•  (Ip'PIpI2)'1  -  2E4|p|'4  +  4m2E4(|pf  |p|4)-‘ 

-  4m2E2E/2(|p'|2|p|4)-1  -  m2o)(|p|2E/)"1  -  frpEflpl’E')-*, 

I*x3  =  4m2E3|p|"4  -  8m4E3(|p'p|p|4)-‘  +  8maE*E/1(|pf|p|4)r> 

+  4m4E(|p'|2|p|2)-1  +  2m2E2(|ppa>)'1  -  2m2E2(|ppa))-1 

-  4m2EE/2(|p'|2|p|2)'1, 


I*,4  -  -2m2E'2ar2  -  2m4E2|p|-4  +  4m6E2(|p'|2|p|4)->  +  4m2E2E'2(|p'p|p|2)-‘ 
-  4m2E/2E4(|J/|2|p|4)-  +  4m4E(|ppai)"1  +  2|^|2m2E2(|p|2a)2)-. 


(691) 


(692) 


(693) 


(694) 


(695) 
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Rewriting  equation  (692), 


!»,.  =  (E|p'|2E/  -  m2co2  -  If/pEoj  +  E'2co2)(frf  |p|2E')- 


or 


U,  =  [ E|p/|2(E/  -  oj)  +  oj2(E/2  -  m2)](|p'|2|p|2E0"’- 
Using  equations  (130)  and  (209)  in  equation  (697),  then 

I»,.  =  EfP|-2. 


Rewriting  equation  (693), 

=  -(|p|2  -  co2)(2|f|2)-  +  I, 

where 


I»,2.  =  -  m2|p|-2  -  4E2m2|p/|'2|p|"2  -  2E4|p|-4  +  4m2E4|^|-2 |p|“4 

-  4m2E2E/2|J/|-2|p|-’  -  m2oj|p|"2E/"‘  -  |£f  E|p|-JE'-. 

Using  equations  (130)  and  (209)  in  the  last  two  terms  of  equation  (700),  equation  (209) 
term,  and  equation  (116)  in  the  fourth  and  fifth  terms,  then 

L«2.  =  -Vi  -  m2|p|"2  -  4E2m2|p/|'2|p|"2  -  2(m2  +  |p|2)2|p|-4  +  4m2E2 
■  (m2  +  |p|2)|H-2|p|-4  -  4m2E2(|p,|2  +  m2)|^|-2|p|- 
-  EE'|p|-2  +  m2|p|-2 


or 


h«2i  =  -Vi  -  4E2m2|p/|'2|p|"2  -  2m4|p|  4  -  4m2|p|'2  -  2  +  4m4E2|p/|‘2|p|'4 
+  4m2E2|p/|"2|p|"2  -  4m2E2fP|-4  -  4m4E2|V|'2|pT  -  EE'IpI’2 


(696) 

(697) 

(698) 

(699) 

(700) 

in  the  sixth 

(701) 

(702) 
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or 


Ux2i - (5/2)(l  +  rn2|p|  2)  -  EE7|p|  2  -  (3/2)m2|p|  2  -  4m2E2|p|  4  -  2m4|p|  4.  (703) 

Then  using  equation  (116)  in  the  first  term  of  equation  (703), 

U21  =  -(5/2)E2|p|-2  -  EE'|p|-2  -  (3/2)m2|p|-2  -  4m2E2|p|-4  -  2m4|p|-4  (704) 

or 

I**2i  =  -(5E2  +  2E7E  +  3m2)(2|p|2)~l  -  2m2(2E2  +  m2)|p|'4.  (705) 

Substituting  equation  (705)  into  equation  (699),  then 

1**2  =  -(Vi)(|p|2  -  o)2)|T|-2  -  (2|p|2)"1(5E2  +  2EE7  +  3m2)  -  2m2|p|-4(2E2  +  m2).  (706) 

Using  equation  (209)  in  the  second  and  third  terms  of  equation  (694),  then 

I**3  =  4m2E3|p|~4  -  8m2(E/2  -  |p'|2  )E3(|pf  IpIT1  +  8m2E3(|pf  +  m2) 

(fpflpl4)"1  +  4m4E(|p|2|p|2)-1  -  4m2EE/2(|p/|2|p|2)*1  (707) 

or 

lbK3  -  20m2E3|p|-4  -  8m2(E/2  -  m2)E3|pt2|p|'4  +  4m2Efpt2|p|'W  -  E72).  (708) 

Using  equation  (209)  in  the  second  and  third  terms,  equation  (708)  becomes 

lbK3  -  20m2E3|p|-4  -  8m2E3|p|"4  -  4m2E|p|'2  (709) 

or 

U.3  =  4m2E|p|'4[2E2  +  (E2  -  |p|2)].  (710) 
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Then  using  equation  (116)  in  equation  (710), 


I**,  =  4m2E|p|'4(2E2  +  m2). 


(711) 


Next,  rewriting  equation  (695), 


i»«4  =  (co2rPi4i^i2)-iFl.4 


(712) 


where 


Filt4  =  -2m2E/2|p|4|p/|2  -  2m4E2co2|p/|2  +  4m6E2co2  +  4m2E2E/2|p|2co2  -  4m2co2E/2E4 

+  4m4co|pfE|p|2  +  2|p'|4m2E2|p|2.  (713) 


Using  equation  (116)  in  the  fifth  term  of  equation  (713),  then 


F»«4  =  -2m2E/2|p|4|p/|2  -  2m4E2co2|p'|2  +  4m6E2co2  +  4m2E2E,2fP|2co2 

-  4m2co2E/2E2(|p|2  +  m2)  +  4m4co|^|2E|p|2  +  2|p'|4m2E2|p|2  (714) 


or 


Ft>4  =  -2m2E/2|p|4|p'|2  -  2m4E2co2|p/|2  +  4m6E2co2  -  4m4co2E/2E2 

+  4m4coE|pf  |p|2  +  2|p'|4m2E2|p|2.  (715) 

Using  equation  (116)  in  the  second  term,  equation  (209)  in  the  fourth  and  sixth  terms,  and  equation 
(130)  in  the  fifth  term,  equation  (715)  becomes 

FbK4  =  -2m2E/2|p|4|p/|2  -  2m4(|p|2  +  m2)co2|^|2  +  4m6E2co2  -  4m4co2 

•  (m2  +  |pf  )E2  +  4m4(E  -  E/)|p'|2H|P|2  +  2|p'|2(E/2  -  m2)m2E2|p|2,  (716) 

or  simplifying  and  using  equations  (116)  and  (130),  then 
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Ftx4  =  -2m2E/2|p|2(E2  -  m2)|p'|2  -  2m6co2fpf  -  2m4fP|2(E  -  E')2|pf 
-  4m4co2|p'|2E2  +  4m4E2|p/|2|p|2  -  4m4|pfEE'|p|2 


(717) 


+  2m2|p'|2E2E'2|p|2  -  2m4|p'|2E2|p|2. 


Equation  (717)  reduces  to 


F.x4  =  -2m6cu2|p'|2  -  4m4cu2|p'|2E2.  (718) 

Substituting  equation  (718)  into  equation  (712),  then 

I.x4  =  -2m4|p|-4(m2  +  2E2).  (719) 

One  next  recalls,  substituting  equation  (691)  into  equation  (676),  that 

1.3  +  1.4  +  1.5  +  U«  =  +  I.«2x"2  +  W3  +  I.«4)C-4.  (720) 

Then  substituting  equation  (720)  into  equation  (578), 

I.  -  I.,  +  I..  +  L,.x-  +  I.x2x"2  +  I.,3x"3  +  I.«4x"4.  (721) 

Next  substituting  equations  (606),  (608),  (698),  (706),  (711),  and  (719)  into  equation  (721),  then 

I.  =  (4|p'|  |T|)-  in  [(|T|  +  |p'|)(|t|  -  |p'|)-][4m2x-2 

-  6 cox"1  -  2co(|p|2  -  co2)x"'|f|"2]  -  x"'|p'|"1  in[(E'  +  |p'|)(E/  -  |p,|)_l] 

+  E'x-IpI"2  +  x"2[-(|p|2  -  co2)(2|T|2)-  -  (5E2  +  2 EE'  +  3m2)(2|p|2)- 

-  2m2(2E2  +  m2)|p|'4]  +  4m2E(2E2  +  m2)x"3|p|"4 

-  2m4x"4|p|"4(m2  +  2E2).  (722) 
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Collecting  terms  with  the  common  factor  (2E2  +  m2),  equation  (722)  can  be  rewritten  as 

I„  =  (2E2  +  m2)|p|'2x'4Gi,  -  (5E2  +  2EE2  +  3m2)(2|p|2jc2)-‘  -  (|p|2  -  co2)(2|T |2jc2)- 
+  E/x',|p|"2  +  (4|p/|  |T|)-  in  ((|t|  +  |p'|) 

'  (|T|  -  |p,|)'l][4m2x'2  -  6o)x"‘  -  2co(|p|2  -  co2)x;_1  |T|-2 ] 

-  jc-Ip'I- in  [(E2  +  Ip'IHE2  -  Ip'I)"1]/  (723) 

where 

G6  =  -2m2x:2|p|“2  +  4m2Ex:|p|"2  -  2m4|p|-2.  (724) 

Equation  (724)  can  be  rewritten 

Gfe  =  2m2|p|_2(-x2  4-  2Ex  -  m2).  (725) 

Substituting  equation  (43)  into  equation  (725),  then 

G„  =  2m2|p|-2[-(E  -  |p|  cos  0)2  +  2E(E  -  |p|  cos  6)  -  m2]  (726) 

or 

Gb  =  2m2|p|_2(E2  -  m2  -  |p|2  cos2  0).  (727) 

Using  equations  (116)  and  (213)  in  equation  (727),  then 

Gb  =  2m2  sin2  0.  (728) 
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Substituting  equation  (728)  into  equation  (723),  then 


l„  =  (l/4){8m2  sin2  0(2E2  +  m2)|p|-V4  -  2(5E2  +  2EE2  +  3m2)|p|'V2  -  2(fP|2  -  oj2) 

•  |T|-2*-2  +  4E/|p|-2k-‘  +  |p/|”‘|t|“1 

•  In  [(T  +  |p2|)(|T|  -  |p/|)‘1][4m2x:“2  -  6ojxT‘ 

-  2co(|p|2  -  co2)|T|-V]  -  4|p'|-,jc-1  In  [(E2  +  |^)(E2  -  l^l)'1]}.  (729) 

Finally,  substituting  equation  (423)  into  equation  (227), 

d3o  =  Z2arf(2n)"1|p/|  |p'|”1m2ar1  dco  d2Q  (Ia  +  I6).  (730) 

Therefore,  substituting  equations  (577)  and  (729)  into  equation  (730)  and  using  equations  (160), 
(161),  (334),  (243),  and  (43),  one  obtains 

d3o/d2Q  dco  =  (Z2e6/8n)(|p/|/|p|co){8m2  sin2  6  (2E2  +  m2)/(|p|V) 

-  2(5H2  +  2EE7  +  3m2)/(|p|2x2)  -  2(|p|2  -  co2)/(|T|2x:2) 

+  4E//(|p|2x)  +  (L/|p||p/|)[4Em2  sin2  6  (3com2  -  |p|2E/)/(|p|2x4) 

+  [4E2(E2  +  E/2)  -  2m2(7E2  -  3EF  +  E/2)  +  2m4]/(|J|2x2) 

+  2co(E2  +  EE'  -  m2)/(|p|2x))  +  Mp'||T|)[4m2/x2 

-  6co/k  -  2co(|p|2  -  co2)/(|T|2x)]  -  4£2/(|p/|x:)}.  (731) 


where 


7c  —  E  -  |p|cos  6, 

(732) 

T  =  p  -  k, 

(733) 

L  =  in  [(EE2  -  m2  +  |p||p2|)/(EE2  -  m2  -  |p||p'|)). 

(734) 

e.  =  <n  [(|T|  +  |^|)/(|T|  -  C'DL 

(735) 

£2  =  In  [(E2  +  |p2|)/(E2  -  IJ2!)]. 

(736) 
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Equation  (731)  verifies  the  Gluckstern-Hull  formula,  equation  (1),  in  agreement  with  equation  (4.1) 
of  Gluckstern  and  Hull's  paper.  1  The  differential  cross-section  equation  (731)  is  that  of  electron- 
nucleus  bremsstrahlung  for  photori  emission  into  solid  angle  d2Q  and  frequency  interval  dco/2rr.  The 
angle  d  is  that  between  the  direction  of  the  emitted  photon  and  that  of  the  incident  electron;  e  is  the 
charge  of  the  electron;  -Ze  is  the  charge  of  the  scattering  nucleus;  p  and  p;  are  the  initial  and  scattered 
electron  momenta,  respectively;  E  and  E/  are  the  initial  and  scattered  electron  energies,  respectively; 
m  is  the  mass  of  the  electron;  co/2tt  is  the  frequency  of  the  emitted  photon;  &  is  the  momentum  of  the 
emitted  photon;  and  T  is  the  difference  in  momentum  between  the  incident  electron  and  the  emitted 
photon.  The  Gluckstern-Hull  formula  is  to  be  used  in  calculations  of  continuum  x-ray  spectra 
radiated  by  anisotropic  relativistic  beam  plasma  systems. 


lR.  L.  Gluckstern  and  M.  H,  Hull,  Physical  Review,  90  (1953),  1030. 
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